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The Larkin-Migdal approach to a cold superfluid Fermi liquid is generalized for a non-equilibrium 
system. The Schwinger-Keldysh diagram technique is applied. The developed formalism is appli- 
cable to the pairing in the states with arbitrary angular momenta. We consider the white body 
radiation problem by calculating probabilities of different direct reactions from a piece of a fermion 
superfluid. The closed diagram technique is formulated in terms of the full Green's functions for 
systems with the pairing correlation. The cutting rules are used to classify the diagrams repre- 
senting one-nucleon, two-nucleon, etc. processes in the matter. The important role of multi-piece 
diagrams for the vector-current conservation is demonstrated. In the case of equilibrated systems, 
dealing with dressed Green's functions, we demonstrate correspondence between calculations in the 
Schwinger-Kadanoff-Baym-Keldysh formalism and the ordinary Matsubara technique. As an exam- 
ple we consider neutrino radiation from the neutron pair breaking and formation processes in case 
of a singlet pairing. Necessary correlation effects are included. The in-medium renormalization of 
normal and anomalous vertices is performed. 

PACS numbers: 
Keywords: 



Contents 

Introduction 

A. Historical remarks 

B. White body radiation 

C. Neutrino cooling of neutron stars 



II. Description of non-equilibrium normal 

Fermi liquids 

A. Dyson equations on Schwinger-Keldysh 
contour 

B. NN interaction and pion degrees of freedom 
in nucleon matter 3 

C. Pion softening and pion condensation ^ 

III. Description of non-equilibrium superfluid 
Fermi liquids 

IV. Optical theorem formalism [l^ 

A. Radiation from a piece of non-equilibrium 
matter [l^ 

B. Diagrammatic decomposition in terms of full 

( — h) Green functions 1^ 

1. Fermions with finite width 

2. Quasiparticle approximation for fermions 1^ 

C. Resummation of the two- fermion interaction 
out of equlibrium. Bosonization of the 
interaction Il9 



D 



Physical meaning of multi-piece diagrams 

1. Non-equilibrium systems 

2. Equilibrium systems 
Extension to a superfluid system 
Application to the point-like interactions 

1. Non-equilibrium systems 

2. Equilibrium systems 



with pairing 
B. Larkin-Migdal equations 



IX. Conclusion 

Acknowledgments 
A. Matrix notation 



V. Renormalization of interaction 

A. Fermi-liquid renormalization 

B. Landau-Migdal parameters for the nuclear 
matter ^ 



VI. Equilibrium systems with pairing at T 7^ Q2 

A. Green functions and response of a system 



VII. Current-current correlators for the PBF 

process n ^ n + I'D [3l| 

VIII. Neutrino emissivity in the neutron PBF 

process [sj 



*The paper is dedicated to the memory of A.B. Migdal on occasion 
of his the 100-th anniversary 
t Electronic address: 



■f Electronic address: 



kolomeitsev@fpv.umb . sk| 
D.Voskresensky@gsi.de| 



B. Lepton tensor 

C. Equilibrium Relations 

D. The loop functions 



2 



References 



I. INTRODUCTION 



A. Historical remarks 



The phenomenological theory of normal Fermi hquids 
at zero temperature was proposed by L.D. Landau in 
Refs. [ll, |2|. A.B. Migdal made the very important ob- 
servation that a jump in the particle momentum distri- 
bution at the Fermi momentum corresponds to a pole of 
the fermion Green's function in the normal Fermi liquid 
and is preserved even in the strongly interacting system 
3 . The presence of the pole contribution to the fermion 
Green's function allowed Galitsky and Migdal to develop 
the Green's function formalism for many-body fermionic 
systems, see Ref. Q. These concepts were first elabo- 
rated on example of low-lying particle-hole excitations in 
Fermi liquids. A.B. Migdal was first who applied these 
methods to description of various nuclear phenomena and 
constructed a closed semi-microscopic approach that is 
usually called " Theory of finite Fermi systems" Q . 

A general understanding of the phenomenon of su- 
perconductivity in case of the weak attraction between 
fermions was achieved by J. Bardeen, L.N. Cooper and 
J.R. Schriffer in Ref. Q, see also Ref. Q for detailed ex- 
position. Due to the sharpness of the Fermi surface pro- 
vided by presence of the Migdal's jump one can consider 
the fermions on the Fermi surface as moving in an effec- 
tive two-dimensional momentum space. It follows imme- 
diately that even a weak attraction between two particles 
is sufficient to form a Cooper pair. As soon as the pairing 
phenomenon is established one can follow different routes 
in description of the superconductivity and superfluid- 
ity: In Ref. N.N. Bogolyubov suggested a very conve- 
nient transformation from the particle ijj operators to the 
new operators of effective excitations on top of the back- 
ground of Cooper pairs. This transformation is broadly 
used in the theory of superconductivity. L.P. Gorkov 
developed the Green's function formalism for supercon- 
ducting fermion systems with an electron-phonon inter- 
action [loj . Y. Nambu introduced a matrix formalism 
to the theory of superconducting metals [111 (Green's 
functions formulated in the Nambu-Gorkov space). In 
Ref. p!] G.M. Eliashberg extended the A.B. Migdal's 
theory of the strong electron-phonon interaction in nor- 
mal metals 13| to include the Cooper pairs. This ap- 
proach can be used to describe strong coupling super- 
conductors. Also A.B. Migdal was the first who rose the 
idea about a possibility of the neutron-neutron pairing 
and superfluidity in neutron stars, which where hypo- 
thetical objects that time [14 1. 

The Fermi liquid theory was then generalized by 
A.I. Larkin and A.B. Migdal for the description of 
fermion superfluids at zero temperature [Tsj . Their for- 
mulation is more general than that done in papers by 
Y. Nambu and L.P. Gorkov, since it allows for differ- 



ent interactions in the particle-particle and particle-hole 
channels. A.J. Leggett applied this formalism for the su- 
perfluid "^He at a finite temperature [l^[l3|- J- Schwinger 
in Ref. [18], L.P. Kadanoff and G. Baym in Ref. [l^, and 
L.P. Keldysh in Ref. [20i | developed the non-equilibrium 
diagram technique for the description of non-equilibrium 
Fermi and Bose systems. Even for equilibrium systems at 
T 7^ Schwinger-Kadanoff-Baym-Keldysh approach is in 
many cases more convenient than the standard Matsub- 
ara technique (applicable only for equilibrium systems) 
since it does not involve the Wick rotation and the ob- 
tained results can be continuously transformed to those 
computed in the standard Feynman-diagram technique 
at zero temperature. 

The importance of coherence time effects on the pro- 
duction and absorption of field quanta from the motion 
of source particles in matter has been first discussed by 
L.D. Landau and I.Ya. Pomeranchuk [2l|. In Ref. [12] 
A.B. Migdal developed the complete theoretical frame- 
work for the description of the bremsstrahlung radia- 
tion of ultra-relativistic electrons in the process of mul- 
tiple rescatterings on Coulomb centers. Successful mea- 
surements of such a suppression of the bremsstrahlung 
radiation have been recently carried out at the Stan- 
ford Linear Accelerator Center [2^, see also the review 
in Ref. Now this effect is named the Landau- 

Pomeranchuk-Migdal effect. 

In the framework of his theory of finite Fermi systems 
A.B. Migdal developed the description of the soft pion de- 
gree of freedom in nuclear matter in application to atomic 
nuclei and neutron stars. In vacuum, pions are the light- 
est quanta of the strong interactions between baryons. 
In medium, pionic modes are softened even further due 
to the coupling to nucleon particle-hole modes and can 
be easily excited even at low excitation energies, simi- 
lar to phonons in solids. As an intriguing consequence 
of the pion softening A.B. Migdal suggested a possibility 
of the pion condensation^ at the increase of the baryon 
density, see Refs. [l^ll^, H^. Latter on, in analogy to 
the pion condensation, the ideas of the kaon condensation 
(29| and the charged p meson condensation [13, HH in the 
interiors of neutron stars were explored. Softening of the 
pionic mode at finite temperature [33-134] and at non- 
equilibrium [35l - [37| may manifest in neutron stars 
and heavy ion collisions [l^, [s^ . 

A.B. Migdal rose question on a possibility of existence 
of superdense abnormal nuclei glued by the pion con- 
densate [1^, H^. Also a possibility of nuclei-stars was 
considered in Ref. [4l|. The similar ideas on a possibil- 
ity of quark nuclei, quark stars and hybrid stars [H, HI] 
are continued to be extensively explored nowadays, see 

MM- 



^ Independently pion condensation was also suggested by D.J. 
Scalapino and R. Sawyer |26|| . 
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B. White body radiation 

Below we consider the white-body radiation from a 
piece of a superfluid fermion matter. To be specific we 
focus on the neutrino radiation from a piece of super- 
fiuid nucleon matter. Standard Feynman technique of 
summation of squared matrix elements of reactions fails 
to calculate reaction rates in the medium, since in general 
case there are no asymptotic states for source particles in 
matter. Indeed, source particles continue to collide be- 
fore and after radiation of a quantum. This gives rise to 
finite imaginary parts of the self-energy functions (par- 
ticle widths). If one naively replaced the summation of 
all perturbative Feynman diagrams (with free Green's 
functions) by the summation of corresponding diagrams 
with dressed Green's functions, it would lead to a dou- 
ble counting due to multiple repetitions of some pro- 
cesses (for an extensive discussion of how one can treat 
this defect see [H, I46l - l48j ). This calls for a formalism 
dealing with closed diagrams (integrated over all possi- 
ble in-medium particle states) with full non-equilibrium 
Green's functions. Such a general formalism was devel- 
oped in Ref. It treats on equal footing one- fermion 
and multi-fermion processes as well as resonance reac- 
tion contributions of the boson origin, such as processes 
with participation of zero sounds and reactions on the 
boson condensates. Decomposition of diagrams is done 
in terms of the full G ^ Green's functions (Wigner densi- 
ties). Each diagram in the series with full Green's func- 
tions is free from the infrared divergencies. In such a 
way one generalizes Landau-Pomeranchuk-Migdal treat- 
ment of the multiple scattering on external centers to the 
treatment of the multiple scattering in matter. Both, the 
correct quasi-particle and quasi-classical limits are recov- 
ered. 

The formulation of the radiation problem in terms of 
closed diagrams calculated within the non-equilibrium 
Green's function in quasi-particle approximation was per- 
formed in Ref. |35|] This approach was called the "opti- 
cal theorem formalism". In Refs. [s^li^ it was demon- 
strated that the standard calculations of reaction rates 
via integration of squared reaction matrix elements and 
the results of the optical theorem formalism match ex- 
actly, provided conditions for the quasi-particle approx- 
imation for fermions are fulfilled. Formally the match- 
ing is done by cutting the closed diagrams. In gen- 
eral case considered in Ref. [47,] going beyond the quasi- 
particle approximation, the series of closed diagrams is 
constructed with respect to the number of the G 
Green's functions. For low temperatures each G ^ line 
brings extra (T/ep)^ factor in the production rate of the 
radiating quanta, ep is the Fermi energy. In Ref. 0| the 
relations between reaction rates at finite widths and the 
quasi-particle rates were found. 

All real calculations of fermion superfluids were per- 
formed within quasi-particle approximation for fermions 
(when fermion width is much less than all other relevant 
energetic scales). Below we focus on the Larkin-Migdal 



approach to the cold fermion superfluids and formulate it 
in terms of the Schwinger-Kadanoff-Baym-Keldysh tech- 
nique to describe fermion superfluids in equilibrium at 
T 7^ and out of equilibrium. 



C. Neutrino cooling of neutron stars 

Physics of neutron star cooling is based on a number 
of ingredients, among which the neutrino emissivity of 
the high-density hadronic matter in the star core is the 
important one. 

After the flrst tens of seconds (at most hours) , the typ- 
ical temperature of a neutron star decreases below the so- 
called neutrino-opacity temperature Topac ^ (1~2) MeV. 
At these conditions neutrinos and anti-neutrinos can be 
radiated directly from the star interiors without subse- 
quent rescattering, since their mean-free path is much 
longer than the star radius [i^. Hence, the star can be 
considered as a piece of a warm "white" body for neu- 
trinos. Typical averaged lepton energy (~ several T) is 
much larger than the nucleon particle width F^v ~ T^/ep- 
Therefore, the nucleons can be treated within the quasi- 
particle approximation. This observation simplifles con- 
sideration essentially. One usually follows an intuitive 
way for the separation of the processes according to their 
phase spaces. The one-nucleon processes (if they are not 
forbidden by the energy-momentum conservation) have 
the largest emissivity, €i, oc for non-superfluid sys- 
tems, then two-nucleon processes come into play, ei, oc 
T^, and so on. In the optical theorem formalism one- 
nucleon processes are determined by the self-energy E 
of virtual W and Z bosons expanded in the series with 
respect to the number, N, of G ^G^ loops with full 
" + and " " vertices. The = 1 diagrams corre- 
spond to one-nucleon processes, the N = 2 diagrams to 
the two-nucleon processes, etc. 

In the so-called "standard scenario " of the neutron star 
cooling, the processes were calculated without taking into 
account in-medium effects. It was argued that the most 
important channel at temperatures up to T ~ 10^-10^ K 
is the modified Urea (MU) process nn — )■ npev. First 
estimates of the MU emissivity were done in [4^, [50| . 
In Ref. [Sll, [S^ B. Friman and O.V. Maxwell recalcu- 
lated the emissivity of this process in the model, where 
the nucleon-nucleon (NN) interaction was approximated 
by a free one-pion exchange (FOPE). The expression of 
the neutrino emissivity obtained by them was used in 
various computer simulations, e.g., in Refs. (53l - [55| . Be- 
sides the MU process, the "standard scenario" includes 
also the processes of the nucleon (neutron and proton) 
bremsstrahlung (NB) nn ^ nnvv and np ^ npvD, 
which contributions to the emissivityis smaller than 
those of the MU processes, see Refs. [5l|, HI]. The den- 
sity dependence of the reaction rates calculated with the 
FOPE is rather weak and the neutrino radiation from a 
neutron star depends very weakly on the star mass. 

There exists another class of so-called "exotic" pro- 
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cesses, which occur only if some special condition is ful- 
filled, i.e. when the nucleon density exceeds some crit- 
ical values. These are the direct Urea (DU) processes 
on nuclcons (e.g., n — >■ peD) and hyperons 57], pion 
Urea reactions on a pion condensate ^46^ 5 8,. 591. kaon 

p-Urca 
DU pro- 
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Urea processes on a kaon condensate 
processes on charged p-condensates [30|, |31| 
cesses on quarks |62| . DU processes on fermion conden- 
sates [1^. The values of critical densities are different 
for various processes and are model dependent. For ex- 
ample, some relativistic mean-field models produce the 
critical density of the DU reaction, n^^ as low as the 
nuclear saturation density no ~ 0.16 fm~'^. However, the 
realistic, microscopically-based Urbana-Argone equation 
of state [i^l yields n^^ ~ 5no. The simulations of the 
neutron star cooling history in Refs. [65l . [66j have shown 
that the occurrence of the DU processes in the neutron 
star with masses M < 1.5Mq would lead to problems 
with the explanation of soft X-ray data. The constraint 
on the equation of state of the dense nuclear matter, re- 
quiring a sufficiently high value of nj?^, was proposed in 
Ref. and explored in details in Ref. [67| . 

It was shown in Refs. [H [H [M El IH Isl that 
the neutrino emission from dense hadronic component 
in neutron stars is subject of strong modifications due to 
collective effects in the nuclear matter. Many new reac- 
tion channels open up in medium in c omp arison to the 
vacuum. In Refs. [11, [H [M El III [11 the nucleon- 
nucleon interaction was considered within the Landau- 
Migdal approach to Fermi liquids. The softening of the 
in-medium one-pion exchange (MOPE) mode and other 
medium polarization effects, like nucleon-nucleon corre- 
lations in the vertices, renormalization of the local part of 
NN interaction due to loop effects, as well as a possibil- 
ity of the neutrino emission from intermediate reaction 
states and DU-like reactions involving zero sounds and 
boson condensates were incorporated. It was demon- 
strated in Refs. [1^, [H, El 113: that for n > no the 
neutrino emissivity is mainly determined by the medium 
modified Urea (MMU) process, in which the neutrino is 
radiated from the intermediate reaction states. This fact 
changes significantly the absolute value and the density 
dependence of the nn — npeP process rate. The lat- 
ter becomes very strong. Therefore, for neutron stars 
with larger masses the resulting emissivity of the MMU 
process proves to be substantially higher than the cor- 
responding value (MU) calculated in the FOPE model 



of Ref. ml. For n > 



no, 



the medium-modified nucleon 



bremsstrahlung (MNB) processes yield a smaller contri- 
bution than MMU ones since the former does not in- 
clude the neutrino radiation from intermediate states, 
e.g. from intermediate pion. However, the MNB pro- 
cesses are more efficient than the NB ones for such den- 
sities. 

Oppositely, for n < no the in-medium effects can 
moderately suppress the two-nucleon reaction rates com- 
pared with those given by the FOPE model [69| - [7l| . The 
pion softening effect disappears at n < 0.5-0.7 ng, see 



Refs. [28|, |38|, but the nucleon-nucleon short-range re- 
pulsion effect remains. Inclusion of the nucleon-nulceon 
correlations without the pion softening [t^I yields a sup- 
pression effect also at n > np. Obviously, this effect also 
follows from general consideration in Refs. (28|,i46,], if one 
artificially suppresses the pion softening effect. 

After the seminal work of A.B. Migdal li*], various as- 
pects of the nucleon superfluidity in neutron stars were 
studied in the literature: The presence of a nucleon su- 
perfluid interacting with the normal component is needed 
for explanation of glitches in pulsar periods and neutron 
star quakes Explanation of pulsar cooling curves 

also requires an inclusion of superfluid phases j66j . Sev- 
eral superfluid phases are found possible. Phase transi- 
tions between different phases may take place fl^ . 

It is commonly accepted that most important are the 
superfluid phases with the spin-singlet pairing of neu- 
trons and protons, in the 15*0 state, and the spin-triplet 
pairing of neutrons in the 3P2 states. The latter is be- 
lieved to occur in neutron star interior at n > ng in the 
state with mj = 0, where mj is the projection of the 
total pair momentum onto a quantization axis. In case 
\mj\ = 2 the exponential suppression of the specific heat 
and the neutrino emissivity is replaced by a power-law 
suppression since the gap vanishes at the poles of the 
Fermi sphere. This possibility was mentioned for the 
first time in Ref. the corresponding reaction rates 
were calculated in Ref. [t^. However, a mechanism to 
realize this interesting possibility in neutron star cooling 
was not elaborated yet. 

Many papers are devoted to the calculation of pair- 
ing gaps within different approaches [76l-l87|. The ob- 
tained results can be essentially different depending on 
a model for the nucleon-nucleon interaction and a cal- 
culation scheme. The predictions of the neutron 3P2 
gaps are especially uncertain, e^., compare Refs. [13] 
and 85 1 . For review see Ref. [83] and references there 
in. Ref. [s^ argues that 3P2 gap should be strongly 
suppressed whereas Ref. [8^1 argues for its strong en- 
hancement. Reference [ill calculated cooling curves us- 
ing both these assumptions and concluded that the cool- 
ing history is naturally explained within assumption on 
the suppressed 3P2 gap. Recently Ref. [11] studied the 
new data on the cooling of Gas A object. Their con- 
clusion is in favor of a suppressed 3P2 gap. At tempera- 
tures below the critical temperatures of the neutron, Tc„, 
and proton, Tcp, pairing, the reaction rates, considered 
above, are suppressed because of a decrease of the avail- 
able phase space. Initially, the suppression effects were 
included simply by multiplying the rate of a two-nucleon 
process by the factors e"^^/-^ [11]. Later, the phase- 
space suppression factors (so called R factors) have been 
treated more accurately in Ref. [ooj . 

In nucleon superfluids, there exist new neutrino- 
production mechanisms, which are forbidden for T > Tr . 
These are the processes, suggested in Ref. [H, [H [9l| . 
in which the creation of a neutrino-anti-neutrino pair is 
associated with the breaking and formation of a Cooper 
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pair - the so-called nucleon pair breaking and forma- 
tion (PBF) processes. The emissivities of the nucleon 
PBF processes are suppressed at T < Tc by the same 
factor - exp(-2A/r) as for the MU, NB, MMU and 
MNB processes. However, in comparison to the all latter 
processes, the nucleon PBF processes have the large one- 
nucleon phase-space volume [s^ H^. The existence of 
this new cooling mechanism demonstrates that influence 
of the nucleon pairing on the neutrino production rates 
cannot be reduced just to an introduction of a simple 
phase-space suppression factor. 

Early works [H, HI, IH, III, which studied the 

PBF processes, did not care about the conservation of 
the weak vector current. The latter is fulfilled only if the 
in-medium renormalization of weak vertex functions is 
performed in accord with the renormalization of Green's 
functions. This problem was tackled in Refs. [93l - l99| . 
Reference [o^ argued that the emissivity of the ISq PBF 
processes should be dramatically suppressed as oc Vp, 
where vf is the Fermi velocity of non-relativistic nucle- 
ons, provided the vector current conservation constraint 
is taken into account. The consistent calculation of the 
PBF emissivity induced by the vector and axial-vector 
currents was performed in Ref. [13, HI] within the Larkin- 
Migdal-Leggett Fermi-liquid approach. The latter takes 
properly into account correlation effects in both particle- 
particle and particle-hole channels. It was demonstrated 
that the neutrino emissivity is actually controlled by the 
axial- vector current and is suppressed only by the factor 
oc Wp, rather than oc Vp. Both neutron PBF and proton 
PBF processes yield contributions of the same order of 
magnitude provided strong and electromagnetic renor- 
malizations of the proton weak vertices llOOl Il0l| 
are included. In Ref. [9^ one argues that for the 3P2 
neutron pairing the vector current conservation changes 
moderately the result obt ained without its inclusion. As 
pointed out in Ref. jl02l |. the suppression of the PBF 
processes at low densities might be served as a possible 
explanation of the superburst ignition. 

As we have mentioned, the convenient Nambu-Gorkov 
formalism developed for the description of metallic su- 
perconductors, cf. Refs. d, [13, [U, does not distin- 
guish interactions in particle-particle and particle-hole 
channels. These interactions can be, however, essen- 
tially different in a strongly interacting system, like in 
the nuclear matter and in the liquid He'^. The adequate 
methods for Fermi liquids with pairing were developed 
for zero temperature by A.I. Larkin and A.B Migdal in 
Ref. [1^ (see also @) and for a finite temperature by 
A.J. Leggett in Ref. [H, |l3|. The problem of calcu- 
lation of a response function of a Fermi system to an 
external interaction becomes tractable at cost of intro- 
duction of a set of Landau-Migdal parameters for quasi- 
particle interactions. Parameters can be either evaluated 
microscopically or extracted from analysis of experimen- 
tal data, see Ref. @. The technical difference of the 
Larkin-Migdal and Leggett approaches is that the for- 
mer approach works out equations for full in-medium ver- 



tices, whereas the latter one calculates directly a response 
function. The former approach was aimed at the study 
of transitions in nuclei, and the latter on the analyzes of 
collective modes in superfluid Fermi liquid. The principal 
equivalence of both approaches was emphasized already 
by A.J. Leggett in Refs. [H, Reference demon- 
strates how one may use both approaches in calculations 
of the PBF rates. 

Reference jl03j | was the first one, in which the most 
important in-medium effects were incorporated in the 
numerical code for neutron star cooling. Among them 
neutron PBF and proton PBF processes were treated as 
equally important. The PBF processes (but with free 
vertices) were i ncorp orated also in the "standar d" co ol- 
ing scenario [13, Il04j that led the authors of Ref. [l05| to 
the suggestion of the minimal cooling paradigm. Detailed 
simulations of different medium effects have been done in 
[65lf66| . In contrast to the minimal cooling paradigm, the 
medium modifications of all reaction rates lead to their 
pronounced density dependence. For the PBF processes 
it is mainly due to the dependence of pairing gaps and 
nucleon-nucleon correlation factors on the density. For 
MMU processes the reaction matrix elements are strongly 
density dependent due to the softening of the exchanged 
pion and the dependence of nucleon-nucleon correlation 
factors on the density. It establishes the strong link 
between the cooling beh avior of a neutron star and its 
mass [M m, m, [el m [103. The density dependence 
of the reaction rates provides a smooth transition from 
"standard" to "non-standard" cooling for the increas- 
ing star-center density, i.e., for increasing the star mass. 
Thus, the inclusion of the most important in-medium ef- 
fects within the "nuclear medium cooling scenario" en- 
ables us to describe appropriately both high and low sur- 
face temperatures obtained from analyzes of soft X-ray 
pulsar data. The mentioned above moderate suppres- 
sion of the PBF emissivity (oc Vp) at 15o pairing should 
not significantly affect general conclusions on the neutron 
star cooling history done in previous works where it was 
not incorporated. 

The paper is organized as follows. In Section |lT] we 
formulate description of normal Fermi liquids at non- 
equilibrium. Softening of pionic degrees of freedom is 
taken into account. In Section IIIII we perform general- 
izations to the fermion superfluids at non-equilibrium. 
The Larkin-Migdal equations are formulated on the 
Schwinger-Keldysh contour. A possibility of the pairing 
in an arbitrary momentum state is considered. In Sec- 
tion [IV] we introduce optical theorem formalism for nor- 
mal and superfiuid fermion systems out of equilibrium. 
Cutting rules for closed diagrams expanded in series with 
respect to the number of G full Green's functions are 
formulated. Important role of multi-piece diagrams is 
shown. Fermi liquid renormalizations are performed in 
Section |Vl Equilibrium T ^ systems with pairing are 
considered in Section IVIl In Section IVIIl as an exam- 
ple, we find the current-current correlator and in Section 
IVIIIl the neutrino emissivity from the PBF processes on 
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to tj t» 

FIG. 1: Closed real-time contour with two external points x, y 
on the contour 

neutrons paired in 15*0 state. Technical details are given 
in Appendices. 

II. DESCRIPTION OF NON-EQUILIBRIUM 
NORMAL FERMI LIQUIDS 

A. Dyson equations on Schwinger-Keldysh contour 

The non-equilibrium theory can be entirely formulated 
on the closed real-time Schwinger-Keldysh contour (see 
Fig. [T|) with the time argument running from tg to oo 
along time-ordered branch and back to along anti-time- 
ordered branch 20]. We assume the reader to be familiar 
with this real time formulation of the non-equilib rium 
theory. Details can be found in Refs. [106 108, lll|,[iil 
and in Appendix A. 

In absence of pairing one deals only with the "normal" 
contour Green's function defined as the expectation value 
of contour-ordered products of operators 

i[Gn{x,y)U = • " 

- <N\Tc^a{x)^l{y)\N> . (1) 

The time-ordering, 7c, goes here according to the time 
parameter running along the time contour C. The av- 
eraging < iV| . . . > is taken over the A^-particle non- 
equilibrium state; a, b are spin indices. In the absence of 
the spin-orbital interaction the Green's function is diag- 
onal in the spin space 

[Gn{x,y)]ab = Gnix,y)Sab, 

where Sab is the Kronecker symbol. 

Because of the two contour branches one actually 
deals with four Green's functions unified in matrices 
(so called Schwinger-Keldysh space). Further explana- 
tions are given in Appendix A. There one can also find 
the helpful relations between Green's functions and self- 
energies. 

The typical interaction time Tint for the change of the 
higher-order correlation functions is usually much less 
than the typical relaxation time Tjei, which determines 
the kinetic stage of the system evolution. Describing the 
system at times t — to ^ Tint , one can neglect initial cor- 
relations. This is in accordance with the Bogolyubov- 
Klimontovich pr inciple o f the weakening of higher or- 
der correlations |l09l Ill0l | . The coarse-graining leads to 
time-irreversibility. Alternatively, one also could suppose 
that the initial state is uncorrelated, like an in equilib- 
rium ideal gas. This corresponds to a situation when 



an inform ation loss occurs right from the beginning, cf. 
Refs. [lol. Illlj . Assuming that we describe the system 
for t 3> Tint we may use the Wick decomposition that 
leads to the Dyson equation formulated on the real time 
contour 

Go^Gn{x,y) ^ Sc{x,y) J dzJ:n{x,z)Gn{z,y), (2) 
where 

Go.'=(^^d, + ^di^St (3) 

in non-relativistic kinematics that we use in this work. 
Here 6c is (5-function on the contour, Go is the free 
Green's function (thin line) 

G^,'G\x,y)^Scix,y), (4) 

and the self-energy I]„ is a functional of the Green's func- 
tions. Being formulated with the standard diagrammatic 
rules the Dyson equation reads 



(5) 



Gn = Go + Go S„ Gn . (6) 

The sign stands here for the contour coordinate fold- 
ing, i.e. the integration, dz, over the spatial coordinate 
and the time coordinate integrated along the Schwinger- 
Keldysh time contour, cf. Eq. (jA2p . Alternatively it 
can be represented as the usual four-dimensional coor- 
dinate integration, if all quantities are treated as ma- 
trices in the Schwinger-Keldysh space, see Eq. (IA7|) in 
Appendix [K\ Thus, Eq. (jG)) is the symbolic equation for 
four (G*-', i,j = — , +) Green's functions. The spin- index 
contractions go in the standard manner: in the direction 
opposite to the fermion arrows. Th roug hout this paper 
we shall use notations of Refs. 0, Ill2j . in which 
and differ by the si gn fro m the corresponding quan- 
tities introduced in Ref. [107[ . In these notations Eq. ^ 
is rewritten as 

Gjf" = Gl' + Gf Ej;' &^ . (7) 

The equation for the retarded Green's function decouples 
from other equations 

Gn = Go + Go G^i , (8) 

and reads in diagrams as Eq. ([5]) (above, this equation 
was formulated for the contour or matrix quantities). 
Here Y.^ — + S ^ is the retarded self-energy, see 
(IA6|) . Similar equation exists for the advanced Green 
function G^. 

We also need to define two-particle Green's function 

Kn{x,y]x' ,y') 

= < N\Tc^-{x')^''{y')^i{x)^{y)\N > 
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and the two-particle interaction amplitude as a con- 
tracted part of K. The system of equations for the non- 
equihbrium t wo-p article Green's functions was studied, 
e.g., in Refs. [lT3l . fTT^ . 

As has been shown by Kadanoff and Baym, in case of 
smooth time-space changes of the system the quasiclassi- 
cal approximation can be apphed to the non-equihbrium 
Dyson equations. Using the first-order gradient expan- 
sion for the quantities in the Wigner representation one 
obtains Kadanoff-Baym kinetic equation for the G 
Green's function, which generalizes the standard Boltz- 
mann kinetic equation for quasi-free particle and the Lan- 
dau kinetic equation for quasi-particle to the case of par- 
ticles with finite mass-widths. This generalized kinetic 
equation is supplemented by the equation for the re- 
tarded Green's function which is algebraic equatio n up 
to second gradients. Self-consistent approximations [ll5j | 
to this ki netic scheme were developed only recently, see 
[l07 L[l08l . [ll6i] and references therein. 



B. A^A'^ interaction and pion degrees of freedom in 
nucleon matter 



Consider the nucleon-nucleon interaction formulated 
within the Fcrmi-liquid approach with the explicit in- 
corporation of the in-medium pion exchange. For zero 
temperature it was done by A.B. Migdal in Ref. [27| . 
and, then, the approach was generalized for finite tem- 
peratures and non-equilibrium systems in Refs. (28l . |32| - 
113, [H, nil . At excitation energies of our interest (e* <C 
ep, e* ~ T in equilibrium) nucleons are only slightly ex- 
cited above their Fermi surfaces and all processes occur 
in a narrow vicinity of the Fermi energy ep. Within 
this approach the long-range processes are treated ex- 
plicitly, whereas short-range processes are described by 
local quantities approximated by phenomenological, so- 
called Landau-Migdal, parameters. At low excitation en- 
ergies the A^A^ interaction amplitude is presented on the 
Schwinger-Keldysh contour (or in matrix notation) as fol- 
lows 



X 



where 




(9) 



(10) 



The solid line stands for a nucleon, the double- line stands 
for a A isobar. Although the mass difference between the 
A and A^, m/^ — mN — 2.1m,r ^ ep ("t-tt is the pion mass) 
the delta-nucleon hole term is numerically rather large, 
since the ttA^A coupling constant is twice larger than the 
ttNN one, and the A spin-isospin degeneracy factor is 4 
times larger than that for nucleons. The doubly-dashed 
line corresponds to the exchange of the free pion with in- 
clusion of the contributions of the residual s-wave ttNN 



interaction and tttt scattering, i.e. the residual irreducible 
interaction to the nucleon particle-hole and delta-nucleon 
hole insertions. The block in Eq. (|10|) . depicted by the 
empty square, is irreducible with respect to particle-hole, 
delta-nucleon hole and pion states and is, by construc- 
tion, essentially more local than the contributions given 
by explicitly presented graphs. Thus the empty block 
term should be much smoother function of its variables 
than the terms (particle-hole, delta-hole, pion) presented 
explicitly in Eqs. (O and (ITUl) . In principle, the short- 
range interaction should be calculated as function of the 
density, neutron and proton concentrations, energy and 
momentum, temperature, etc. However, instead of doing 
complicated calculations one often reduces it to the set 
of Landau-Migdal parameters, which one extracts from 
analysis of experimental data on atomic nuclei. 

The irreducible part of the interaction involving A iso- 
bar is constructed similarly to (ITU|) 




>■■( 



(11) 



The main part of the A^A interaction is due to the pion 
exchange. Although information on local part of the NA 
interaction is rather scarce, one can conclude [H, Ill7j | 
that the corresponding Landau-Migdal parameters are 
essentially smaller then those for A^A^ interaction. There- 
fore, and also for the sake of simplicity we will, in further, 
neglect the first graph on the right-hand side of Eq. (|lll) . 

The spectrum of the particle excitations is determined 
by the spectral function given by the imaginary part of 
the retarded Green's function {A = Resumma- 
tion of diagrams shown in ([9|) yields the following Dyson 
equation for the retarded pion Green's function 



(12) 

Here E^g^ is the residual retarded pion self-energy that 
includes the contribution of all diagrams, which are not 
presented explicitly in (|12|) . like s-wave nN and tttt scat- 
terings (included by doubly-dashed line in ([rUl) ). The full 
vertex takes into account A^A^ correlations 






(13) 

In some regions of the (w, k) plane the pion spectral 
function A^(a;, k) has sharp peaks along some lines U!i{k) 
which we call the spectral branches. Nearby these lines 
one can use the quasi-particle approximation^ and write 



^ The quasi-particle approximation for the given particle species 
is understood as putting the imaginary part of the retarded self- 
energy in the Green's function to zero. The quasi-particle width 
r = — 2QE^ is then calculated with so-defined Green's functions. 
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FIG. 2: The typical behavior of the vr*'" spectra in symmetric 
nuclear matter and of tt" spectrum in asymmetric matter. 
The chosen typical value of the density n = 2no is supposed to 
be smaller than the critical density for the pion condensation. 

the spectral function as 

27r 5(a; — a;i(fc)) , , 

* I- L.(fc) 

The spectrum of pionic quasi-particles possesses three 
branches for tt^'" in the symmetric nuclear matter {N — 
Z) and for tt*' in asymmetric matter (N ^ Z), e.g. in 
the neutron star matter. The typical spectrum is shown 
in Fig. [2j In the region uj ^ to^ there are two branches: 
the A branch and the pion branch. For uj < TOtt there 
is the spin sound branch (with w — when k ^ 0). 
The hatches show the regions on the (w, k) plane with 
a non-vanishing pion width, calculated within the quasi- 
particle approximation for nucleons and A isobars. In 
the lower hatched region, at cj < kp^.N f^nd k ^ Pf.at, 
there are no quasi-particle branches and the pion width 
cannot be neglected. This is the region of the Landau 
damping in the nucleon particle-hole channel. The pion 
spectral function is enhanced in this region of uj and k 
for n > Hci ^ (0.5-0.7)no. We stress that the pion spec- 
tral function calculated beyond the quasi-particle approx- 
imat ion f or nucleons and A's is much more involved, see 
Ref. [nil- 

To specify the enhancement of the tt*'" spectral den- 
sity for N = Z and of the 7r° for N ^ Z in the Landau 
damping region it is convenient to introduce the function 

^^{k) ^ -[m'^iuj ^ o,k, 

Note that momenta passing through the NN interaction 
in the MU and MMU processes are k ~ PF,n, where pF,n 
is the nucleon Fermi momentum, and for the MNB pro- 
cess k ~ (0.9-1) pF,n (46] . Remarkably, the minimum on 



the function uj'^{k) is realized at the similar momentum 
ko ~ (0.9-1) pF,ji- Thus, the quantity u!{ko) called the 
effective pion gap controls the strength of the A'^A'^ inter- 
action. The NN cross-section is oc 1/w'* for tj^ < mj, 
provided the MOPE dominates, see Eq. below. Note 
that for the asymmetric nucleon matter the pion gap is 
different for tt" and for tt^ since neutral and charged 
channels are characterized by different diagrams permit- 
ted by the charge conservation. 

The pion chemical potentials (/i^+ ^ /^^- ^ 0, = 
0) are determined from equilibrium conditions for the 
reactions involving the corresponding pions. In the neu- 
tron star matter /z^- follows from the condition of the 
chemical quasi-equilibrium with respect to the reactions 
n — pTr~ and n — pev: /z^- = He = £F,n — £F,p, where 
£F.n, EF.p are Fermi energies of the neutron and proton. 
For a small-size systems like atomic nucleus one should 
put = fl^- = fl^o = 0. 

At low pion energies (for 7r*^° for iV = Z and for tt" for 
N Z) the lowest-energy state determining by the pole 
of the pion Green's function is i/Su ~ a;^(A:o) with /3 > 
appeared due to the Landau damping, see Ref. [28l| . Thus 
for tj^(fco) > the pion excitations die out with time 
exponentially oc exp(— w^(fco) t//3). 



C. Pion softening and pion condensation 

For w^(fco) < the pion field grows exponentially with 
time as exp(|w^(/co)| t//3). Thus, the change in the sign 
of cD^ at n = nc,,r marks the critical point of a phase 
transition to a state with a classical pion field (a pion 
condensate). The critical density Uc^tt depends on the 
values of Landau-Migdal parameters, which are badly 
known for asymmetric matter and for densities signifi- 
cantly larger than uq. Nevertheless, some estimations can 
be given. Various experiments have shown that the pion 
condensation does not manifest itself in atomic nuclei 
as a volume effect, see Refs. [13, Hl]^. Different model- 
dependent estimations indicate that Uc^ti ^ (1.5 — 3)no, 
depending on the pion species, the proton-t o- neutron ra- 
tio and the model used, see Refs. pg-ITTl. Variational 
calculations |l2l| yield Uctt — 2no for isotopically sym- 
metric nuclear matter and Uc^tt — 1.3no for tt^ mesons in 
the neutron star matter. 

Typical density behavior of w^(fco) (for 7r^,7r° at iV = 
Z and for 7r° at iV ^ Z) is shown in Fig. [31 At ti < rtci , 
uP{k) has the minimum for fcg = 0, i.e. uj^{kQ = 0) = 
rn^ — ^j^. For such densities the value u;^(pF,n) essentially 
deviates from — tending to + Pf « ~ A '^^ the 
low density limit. 

At the critical point of the pion condensation (n = 



In the surface layer of the nucleus a classical pion condensate 
field may appear due to a coupling of the gradient of the a moan 
field to the pion, see |12C| . 
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FIG. 3: Effective pion gap (for 7r±'" for iV = Z and for tt" for 
N ^ Z) versus the baryon density as calculated in Refs. [H, 
|38| |. Compare with the pion spectrum presented in Fig. [2] 



nc,7r) the value ^^(fco) with artificially neglected tttt fluc- 
tuations changes its sign (dashed line in Fig. [3]) . It sym- 
bolizes the occurrence of a second-order phase transition 
to an inhomogeneous (feg 7^ 0) pion-condensate state. In 
reality, the tttt fluctuations are significant in the vicinity 
of the critical poin t and the phase transition is of the first 
order [sl, HI, Il22l | . Therefore we depict two branches in 
Fig. [3] (solid curves) with positive and negative values of 
w^(fco). Calculations in Ref. [s^l demonstrated that at 
n > ric^TT the free energy of the state with a;^(fco) > and 
without the pion mean field becomes larger than the free 
energy of the state with uP(kQ) < and a finite pion mean 
field. Therefore at n > nc,7r the state with cD^(fco) > is 
metastable and the state with oP{kQ) < and the pion 
mean field <y9^ ^ becomes the ground state. 

The quantity uj'^{ko) demonstrates how much the vir- 
tual (particle-hole) mode with pion quantum numbers is 
softened at the given density. For the symmetric nuclear 
matter at n — no the ratio a — G-^/G'^ ~ 6 for w = 0, 
k = Pf,n- However, this so-called "pion softening" (27j 
does not significantly enhance the NN scattering cross 
section because of the simultaneous essential suppression 
of the ttNN vertex by nucleon-nucleon correlations. In- 
deed, the ratio of the NN cross sections calculated with 
the FOPE and MOPE equals to 



R = 



cr[FOPE] _ -f^{Lo ~ 0, fc ~ PF,N){ml + pI j^Y 
a[MOPE] ~ 



Kpf.n) 



(15) 

where 7 is the vertex dressing factor determined by 
Eq. p^ . 7(^0) — 0.4. For n < uq one has i? < 1, whereas 
already for n = 2no this estimate yields R ~ 10. Thus, 
following Refs. 27, 28] one can evaluate the NN inter- 
action for n > no with the help of the MOPE, i.e.. 



(16) 



Here the bold wavy line relates to the in-medium pion. In 
the soft-pion approximation the same one-pion exchange 
determines also interaction in the particle-particle chan- 



X 



(17) 



Namely, this quantity determines the NN interaction en- 
tering neutrino emissivities of the MMU and MNB pro- 
cesses. 

Often, one considers the softening of the one-pion ex- 
change only, neglecting the same effects for the two, three 
etc. pion exchanges, arguing for their smallness because 
of extra integrations over the intermediate states, see 
Ref. [27} . At zero temperature these effects are numeri- 
cally small. Nevertheless, they result in the change of the 
order of the p hase transition (from the second order to 
the first order) (l22l |. Ignoring this small jump in the pion 
gap one may deal with the MOPE interaction in the one 
(particle-hole or particle-particle) channel. The residual 
interactions are then hidden in the values of Landau- 
Migdal parameters. In case of the non-equilibrium and 
equilibrium T ^ matter these pion fluctuation ef- 
fects contribute significantly to the pion self-energy and 
should be taken into account, see Refs. [UHJ]. Corre- 
sponding pion fluctuation contributions must be then ex- 
tracted from the residual interaction parameters. Their 
self-consistent analysis in both particle-hole and particle- 
particle channels has been performed in Ref. [3J| within 
the Thomas-Fermi approximation (fco ^ 2pF,Ar)- 



III. DESCRIPTION OF NON-EQUILIBRIUM 
SUPERFLUID FERMI LIQUIDS 

Below, dealing with pairing phenomena in non- 
equilibrium systems we assume that deviations from the 
equilibrium are rather small, bearing in mind that non- 
equilibrium effects should not destroy the fermion pair- 
ing. 

In a superfluid Fermi system a condensate of paired 
fermions is formed. It induces non- vanishing amplitudes 
for the transitions of a particle to a hole state and vice 
versa. Thus, it is possible to combine one-particle state 
on top of the iV-particle background and one-hole state 
on top of the background with N + 2 particles. The one- 
particle-one-hole irreducible amplitudes of such transi- 
tions are depicted as blocks 

-z[A«(a;)]ah = -I <N\^a{x)^b{x)\N + 2> 



— b- 

^i[A^^\x)U - -I <N + 2\^i{x)^{x)\N> 



(18) 

Their spin structure can be written in general case as 

A(2) (x) = ta2 ( A^'' (x) ao + ^ {x) a) . (19) 
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Here a = (cri, (T2, ua) with Gi being the Pauh matrices; 
ctq is the unit matrix in the spin space. The amplitudes 

Aq ' are induced by a particle-particle interaction with 
even angular momenta L — 0,2,4,..., and amplitudes 
Aj^ ' , by the interaction with odd angular momenta L = 
1,3,5,..^. 

Since A^^^^^ are functions of the only one contour coor- 
dinate, we have [A(1'2)]-+ = [A(1'2)]+- = and denote 

[^(1,2)]— ^ [A(1^2)p ^ -[A(l'2)]++^ foUo^g f^Oj^ 

Eq. (jA5[) valid for any two-point functions. From the 
definitions (IT51) follows 



(20) 



(a[,^V = a^^\ {K^l^y^-Kf\ (21) 

The transition amplitudes (jl8l) imply that a propagat- 
ing particle can be transformed in flight into a hole and 
vice versa. This is described by a new type of propa- 
gators — anamalous Green's functions — which can be 
defined on the Schwinger-Keldysh contour as 



(A(i)(x))t = (A(2)(x))T 
and Eci. (HH) yields 



^<N\Tc^a{x)^b{y)\N + 2>, 

= <N + 2\Tc^iix)^iy)\N> . 

For these two functions are valid relations 

[F(^'^^\x,y)U = -[F^'^'Hy,x)]ba, 
[F^'''Hx,y)]l = -[F^'-''Hy,x)U. 

We also need the Green's function for a hole 

[iG^'{x,y)]ab = a —6 

=<N\Tc^cAx)% ,{y)\N >, 



(22) 



(23) 
(24) 



(25) 



which is written in terms of the charge conjugated 

fermion operators \l/c = 0^2 (^^) j where T stands for the 
transposition operation in the spin space. The Green's 
functions of a particle and a hole are related as 



G\x, y) = fT2 [G(y, x)]^ = G{y, x) Sab ■ 



(26) 



The Dyson equations for dressed normal and anoma- 
lous Green's functions in case of pairing (Gor'kov equa- 
tions) have the form: 



(27) 
(28) 
(29) 



The thin line is the normal free Green's function of the 
fermion particle and the inverted thin line is the normal 
free Green's function of the hole, E is the full self-energy 
of the particle and E'', of the hole. Explicitly these Dyson 
equations read as 

G==Go + Go0E0G-t-Go0 A^i) F^^) ^ 
F^^') = E'' +G^Q A(2) G , 
/'(I) = Go E F^^) -f Go A(i) G'' . (30) 

In terms of the dressed Green's functions for the system 
without pairing ([6]) the Gor'kov equations can be written 
shortly as [15], 

G~G„ + G„0A(i)0F(2)^ 

f(2) ~ G'^, A(2) G , 

~ G„ A(^) G'' . (31) 

In these equations we dropped the terms containing the 
differences (E - E„) and (E'' - Ejj). The self-energy 
E includes the same set of diagrams as the self-energy 
E„, but constructed now from the Green's functions for 
the superfluid system, G and F, instead of the normal 
Green's function G„. Note that the anomalous Green's 
functions F can enter E only in pairs to preserve the 
number of incoming and outgoing fermion lines in each 
vertex. Hence, the terms containing F's are oc A^. Since 
in the momentum representation G — G,i oc A^ and G — >■ 
Gn for Ipo — cfI/A >> 1, each integral over the energy 
Po in the self-energy is accumulated only in the vicinity 
of the Fermi surface. Thus the neglected terms are small 
as AVe|, cf. also Ref. O. 

Hereafter the thin line in diagrams will stand for the 
full Green's function of the system without pairing. From 
Eq. (PT|) it follows immediately that the Green's function 
G is diagonal in the spin space. Then the Green's func- 
tions F*^^^ and F*^^) have the same spin structure as the 
blocks (UHl): 

F«(a;,y) = {F^'\x,y) ao + Fl^\x,y) <j)ia2 , 
F<^^Hx..y) - i<J2{F^''\x,y)aQ + F^^\x,y)a). (32) 

In further, we will assume that only one type of pairing 
may occur: either in a state with an even angular mo- 
mentum or in a state with an odd angular momentum.^ 

Then wc have either only Ag^'^-* ^ or only A^^'^-* ^ 0, 

(1 2) ^ (1 2) 

and correspondingly, either Fq 7^ or Fj^ 7^ 0, 

but not both simultaneously. In this case the Green's 

function G remains diagonal in the spin space. Then we 

may follow derivations of [15] with a little difference that 



^ It is however not excluded that at some conditions a part of 
fermions can be paired in one state whereas other part, in another 
state. 
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all our expressions are valid on the Schwinger-Keldysh 
contour. 

The graphical equations for A^^^ and A^^^ are 



(33) 



4^ 



(34) 



Here is a two-particle irreducible interaction, that de- 
termines the full in- medium particle-particle ( "pp" ) scat- 
tering amplitude via graphical equation 




^^^Xltt^' ^^^^ 



which in terms of contour foldings becomes 

fp^ = V -lVU]{Gn Gn) □ fpp . 



(36) 



Both V and T are to be understood as formulated on the 
Schwinger-Keldysh contour, 



Lpp 

-i 



dzi I dz2 / dz3 / dziV{x' ,y';zi,Z2) 
c Jc Jc Jc 



xG'„(zi, Z3) G„(Z2, 2:4)Tpp(z3, Z4; X, y). 



(37) 



or being matrices in the Schwinger-Keldysh space. The 
sign □ in (1361) stands for integration over two four- 
dimensional coordinates with the time running over the 
Schwinger-Keldysh contour. The quantities V and T are 
also matrices in the spin space, and we parameterize V 
as 



Vcd,ab = i 



(38) 



= Vo{i(72)dc{i(T2)ab + V"'^[ff°'ia2 )dc{icr2^'^)ab ■ 

With this definition the interactions Vq and Vi corre- 
spond to the scattering of two fermions with the total 
spin zero and one, respectively. The interaction block 
entering Eq. ([55]) is given by 



Vab,cd ■ 



(39) 



Using this definition we can analyze the spin structure of 
Eq. dSl) 



[A^'%^-l[F^%daV,d^ab 

= -i[Gji]cc [A^'^k Gdd □ Vcd.ab ■ (40) 
Substituting here Eqs. (|19|32p and (|38|) and taking into 
account that the full Green's functions G and G„ are 
diagonal in the spin-space, we obtain 



= -I [G':, za2(A[,') + a) G]^, 

^\Vo{i(T2)dc(i^2)ab + V;"''((T"icr2)dc(icr2(T'^)a6]- 

Separating terms with ctq and a Pauli matrices and mak- 
ing use of the relations {10-2)^ — —1 and Tr{a"(T^} = 
2 S'^l^ we find 

a[,'' =2z(Gf,0A[,') 0G)Hyo, 

Af = 2 z (G^ A'-f^^ G) □ 1//" . (41) 



For Aq '' and A^^'' we obtain exactly the same equations 



and aW 

as (|^T|) but with the replacement G^ ^ G. 

An external field V°^*^ acting on a supcrfluid Fermi 
system can induce four different effects determined by the 
four vertex functions related to the creation of particle 
and hole r, anti-particle and anti-hole t'*, two particles 
r'^-', and two holes t'^-*. The vertices can be graphically 
depicted as 



-IT 



(1) 



r(2) 



(42) 



In matrix notations, vertices have three indices, t'*'' 



where the lower index relates to the external dash line. 
The couplings of the external field to the particle and to 
the hole are related as 



r''{x,z,y) = [f(y,z,a;)]'^ 



(43) 



The coupling of an external field to the non-rclativistic 
fermion is described by the 2x2 matrix acting in the 
fermion spin space. Any rank-2 matrix can be decom- 
posed into a unit matrix do and Pauli matrices a. Thus, 
we have 
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f — to ao + a ti , = tQ (7a + a'^ . 



r(2) ^ 



ia2 (tg ctq + O" ^1 



(2)^ 



and similarly for the bare vertices 



To — tofif^O + (^toA , To = to.0^0 + to^l 



The vertices obey the following equations defined on the Schwinger-Keldysh contour: 



(44) 



(45) 



(46) 

Intermediate lines in (I46p are of all possible signs in the Schwinger-Keldysh space. In operator form equations for 
vertices read 

f = fo~ i\GQTQG + f'' + g q 7^(1) ^(2) + p{i) Q ^(2) cj □ {7 ^ (47a) 
G'' f'' G'' + q ^ q pW 4_ ^(2) Q ^(1) G'» + G'' f (2) ^(1)1 □ u ^ (47b) 
G f G'' + F^^^ f(2) /'(I) + G f + f'' G''j □ T> , (47c) 
i G'* f(2) G + Z'f^) ^(1) p{2) ^ ^(2) 07-0G_,_G'j0 7./i0 p{2)] □ ^ ^ (47^) 



f(i) 



f(2) 



f ^ - i 



r 



Here C/ is the particle-hole irreducible interaction, which 
determines the full particle-hole scattering amplitude via 
the equation 



fph = C/-zC/H(GG")Hfph. 
In diagrams this equation has the form 



(48) 




^^XlK^' ^^^^ 



which differs from (1351) by inversion of one of the nucleon 
lines. In the spin space the matrix U is defined as 



Udc,ab = i 



UoSdcSab + ufaia^, 



(50) 



The amplitudes Uo and Ui are determined on the contour 
or they are matrices in the Schwinger-Keldysh space. 



Having at our disposal the compete spin structure of all 
elements we can work out the spin algebra in Eqs. (j47p . 
As an illustration we consider the second and third terms 
in Eq. (H7a| : 
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+ Gao (4'^ ffO + ICT2 iCT2 (F^'^ ^70 + ^ C?) + . . . ] □ (C/o ^<ic<5a6 + uf af, <b) . 

For the spin-scalar to component we obtain 



to - to,o = -2^[- (^(1) t^, F^'^ + F^'^ ^ Ff )) - G (t« F^^' + T/^) f/^^) + 



■ 



(51) 



fl 2) '^(12) 

We dropped here ah terms containing simuhaneously Fq ' and F^^ Green's functions, since we assumed only one 



type of pairing in our system. Using the trace Tr{(T"cr^cr^} — i 2 Eq^^ we obtain the term iea^-yP}'''^"' 1^'"^^^ 



!(2)7 



which vanishes since the vectors F^^) and F*^^^ are colinear and may differ only by a phase, see Eq. (|24 
The equation for the spin-vector vertex ti is more involved 



- G (r/^)'^ f(^) + {ti'^Sf^s + ^e0,st}'^'') f/^'*) +...]□ C/f^ (52) 

We used here the trace TT{a"a^a^a^} — 2 (SapS^s — Sa^Sps + SasSp^) . The different signs of terms in the first bracket 
appear because 172 i?"^ (72 = —a. We elaborate the spin structure of V at hand of the second and fourth terms in 
Eq. B7a : 



ab 



,(2) 



(2)^ 



+ (Fq'^' ao + Fi^'^ a)ia2 (t^ ao + a^ T^) GVo + . . ■] jVo{ia2) dSa2) ab + V^i"''(«a2a'^)dc(a"*a2 )afc) • 
This equation reduces to 



''0 



= -2i [ 
= -2i\ 



+ (f(i) f ) Fo^^^ -f F^'^ F|^0 



(l)^i(2) 



©4 



Fi(')0i7')0G'' + .. 



+ {F','' i?^^'^ F„^^^ - 6,x5sp + 6,p6sx)F^'>' i?^^' F^''^) 

- (f/1)^ {t'^S.p - iep,stn - F^i) f/'^) G" -f . . . ] □ V;"^ 



(53) 
(54) 



Proceeding this way and collecting all terms we rewrite Eqs. (j47p for the scalar vertices as 



to — to 



— f 



0,0 



-2i 



~2i 



GQtoQG- F^^'^ 0^0 F^ 



n(2) 



G0 4^)0Ff 



fW0 4^)0G 



G 4*^^ F^;'' - Fl'> ^1^"^ G 



;(2) 



(2) 



'(1) 



G'' To'' G'' - Fq^') to F^^^ 



F 



(2)^,(1) 



Ff 0to0F(^^ 



f/'^ ii^^^ G 



(1) 



0*0 

h 



QG'' 



□ f/o, 
^0 QF^' 



G"0 42)- 



G'' 4^'^ f/'^ HJ/o, 



(55a) 



(55b) 



4'^ = -f 2 i [G 4'^ - Fq^') 4'^ Fq^') -t- G to fI^^> + fI^^> Qt'^QG'' 

- 4'^ ^1'^ + G Ti f4'^ - f4'^ ti C*] □ Fo , (55c) 
4'^ = +2 i [G'' 4'^ G - Fq*') 4'^ Fq*'^ + Fq^"' to G + G'' 4 F^^^^ 

- Ff ^ 4'^ ^ + ^/'^ fl G - G'' fl'' f/'^] □ Fo . (55d) 
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For the vector vertices we have 



''0,1 



7*ha 
'-0,1 



-(l)a 



2 i [G fi" G + Fq^^' Ti'"^ 



^.(2) 



G0r«'^0Ff 



F^'^ ti^^^*^ G 



(2)/3 , 



- G (t«<5^5 + lep.stT^-') f/')' - f/^)^ {t^^^S^p + lefs.stl'^') G] □ Uf . 



G'^ r/"" G" + Fo^^^ fP Fo^^^ + F^^' ^0^^^-^ G'^ + G'' io^^^'^ F, 



p(2) 



(2) 



r(i)/3 , 



-(2)/3 , 



+ {S^sSxp - S^xSsp + (S^/3(55a)f/^''' Fi 

+ 2 ^ [G r (1)^ G" - f(^) r/^)^ F^^) + G F^'^ - F^'^ t,^^ G 



(1)A 



(56a) 



(56b) 



~ (<5^5<5a^ - S^xSsp + S.,p5sx)F['^'' f/'^* F^^^^ 

+ G (to<5^5 + lep^si^) f/^)' + F^'^^^ (t^JJ^^ - le^^^f^''^) G"] □ F^"^ , (56c) 

r/'^" = + 2 z [G'' r (2)^ G - f(^) r/^)'' f'^^ + f'^^ g - g'^ 1,'-^ f*^) 

- (<5^5<5a;3 " S^xSsp + 5^pSsx)F?^'' t^^^' Ff 

+ i^/'^'' {to5^p + G + G" {t'^^Sps - le^^stl"'') f/'^*] □ F^"^ . (56d) 

Recah that Eqs. (1551) and (|56p are written here simuhaneously for both types of pairing with even and odd angular 

(12). . 2) 

nomenta. For the former case we must retain only the terms with Fq ' , in the latter one, the terms with F^ ' . 
Response of the system to the external field is described by the self-energy 

G ? G + F(i) f F(2) + G f F(2) + f(i) f (2) c 



= -w 

After taking the spin trace we obtain 



(57) 
(58) 



S = -2ito,o □ 



G to G - F(i'^ tj Fo^'^ - G F(l'^ - F(l'^ G 



p(2) 



(1) 



r(2) 



f/^) F/^^ - G F,''> - F,''' t7^^ G 



(2) 



r(i) 



(2) 



(1) 



^(2) 



-2 z to"i □ G ti" G + Fo^'^ ti'*" Fo^"^ ~ G t^^'^" F^^> - F^^^ t^'^>" G 



-(2)a 



G (4^)<S„, + Z6„^,t7^'0 F/^^-^ - F/^^'^ {t^'>S^0 + ze„^,r/^^^) G 



(2)7 



(l)/3 



(59) 



If external field couples to a conserved current the self- 
energy (|58p must support the current conservation obey- 
ing the relations 



d^i:^,{x,y) = d;j:^.{x,y) = 0. (60) 



Following our notations the first Lorentz index, fj,, in 
T,fj_^{x,y) is attached to the right vertex (at the contour 
coordinate x) in diagrams (|58p and the second index, v, is 
attached to the left vertex (at the contour coordinate y) . 
Relations (pD|) can be fulfilled if the full vertex functions 



p2)) satisfy the Ward identity 
G{x,z) - G{z,y) 

= G{x,x') Q d^T^{x\z,y') Q G{y\y) 
+F(i) {x, x') d>^,?l {x',z, y') F(2) (y', y) 
+G{x, x') 9,^r(i) (x', z, 2/') F(2) (2;', y) 
+F(i)(x, x') a,^f f [x', z, y') G(y', y) (61) 

derived first in Ref. 15] and generalized here for coordi- 
nates on the Schwinger-Keldysh contour or for matrices 
in the Schwinger-Keldysh space. 
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IV. OPTICAL THEOREM FORMALISM 

A. Radiation from a piece of non-equilibrium 
matter 

Optical tlicorcm formalism is an efficient tool to cal- 
culate reaction rates including; finite particle widths and 
other in-medium effects [s^ 113, 111]. Assume that we deal 
with a system of a finite size (white body) transparent 
for radiating quanta. To be specific let us consider anti- 
neutrino-lepton (P/) production. By the lepton we mean 
the electron e, muon /i^, or neutrino v. We assume that 
the system is opaque for e and /i^ but transparent for 
V and V. Then it is convenient to express all quantities 
in the Wigner representation doing the transformation 
(1X9)) from coordinates x = X - ^/2, y = X + ^/2lo the 
corresponding Wigner coordinates {X]q) = {t,X;uj,q). 
The probability of the anti-neutrino-lepton production 
can be expressed in terms of the evolution operator S, 



d>V^°i. 



qi - qp) 



d^xdtd^q 

-iT.-+{x,q) = J2 {'^l\SS^\X*){X*\SS\Dl) , 



where dS = S — I, and 



d$. 



(1 - ni)dqfdql 
(27r)6 4w; Up 



(62) 
(63) 

(64) 



is the phase-space volume of an antineutrino with the 
four- momentum qp = {ujp,qp) and a lepton with the 
four- momentum qi — {uji,qi). If / = v, the occupation 
number ni is to be put zero. The summation goes over 
the complete set of all possible intermediate states {^*}. 
In Eq. ([63)) we suppose that electrons or muons can be 
treated in the quasi-particle approximation. Then there 
is no need (although possible) to consider them as inter- 
mediate states. 

Making use of smallness of the weak coupling, we ex- 
pand the evolution operator as 



+ 00 



SS'. 



T {Vwixo) Snucl{xo)}dxo , (65) 



where Vw is the Hamiltonian of the weak interactions 
in the interaction representation, S'nuci is the part of 
the S matrix corresponding to strong nuclear interac- 
tions, and T{...} is the chronological time ordering op- 
erator. After substitution of 5S into (|63)) and averaging 
over a non-equilibrium initial state of the nuclear sys- 
tem, there appear chronologically ordered (G , F ) , 
anti-chronologically ordered {G~^^ , F++) and disordered 
{G^ , and G ^, F ^) exact Green's functions. 

Once the reaction probability is evaluated according to 
Eq. (1621) . the neutrino emissivity in the neutral channel, 
i.e. with Pi' production, is given by 



d>V^°i_ 



d^xdtd-^q 



CO duj d q . 



(66) 



The emissivity in the charged channel, i.e. with vl pro- 
duction (Z = , /i~), is as follows 



dW^°i„ 



d^xdtd^qd^qp 



Wyd^gpd q. 



(67) 



The Lagrangian density for the lepton-nucleon inter- 
actions is 

G 

^ (.?ch,^^ch + Jnout,;.Cut + h.C.) , (68) 

where G = 1.166 x 10~^ GeV~^ is the Fermi coupling 
constant and there are two contributions from charged 
(ch) and neutral (neut) currents; "h.c." stands for her- 
mitian conjugated terms. The lepton currents are 



neut 



^ch = V'e7^(l - 75)^. 



(69) 



The nucleon currents jch and jncut have vector and 
axial- vector components 

Jch.fj. — I Jncnt.fj. — ^ J fj. i 



jT = -T;<iP%J^ - 9A J^)^n{p). (70) 



where the four vectors J^'^ c an b e written for non- 
relativistic nucleons as, cf. Ref. |l23l |. 



J' 



(1, f±2) , J^ = { 



2miv " 



(71) 



and Ci, = 1 — 4sin^0w — 0.08, 6'w is the Weinberg 
angle; qa — 1.26 is the axial-vector coupling constant, 
P — ipoiP), P and p' are momenta of incoming and out- 
going nucleons. Compared to the frequently used ex- 
pression, that includes only leading terms in the non- 
relativistic expansion, e.g., see Ref. |5lJ, we following 
[9^ [95} retain here sub- leading terms oc vp . Although 
in many cases the terms cx vp yield small corrections 
to leading-order results, in some particular cases the 
leading-term contribution may vanish because of sym- 
metry constraints and then sub-leading terms become 
dominant. Such an example will be studied below. The 
bare current vertex involves the bare nucleon mass, see 
Eq. (j71l) . In medium, however, nucleon wave functions 
are normalized to one quasi-particle rather than to one 
free particle, provided the nucleons are treated in the 
quasi-particle approximation. Hence, the bare nucleon 
mass mjv is to be replaced by the in-medium nucleon 
mass TO^. 

The structure of the weak-interaction Lagrangian ()68p 
suggests that we can detach leptonic currents 



(^2 

2 nuc\, fiv / J L J ' 

spin 



(72) 
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and deal with the object determined by the strong inter- 
actions only 

d^^e^'^^iUUX - e/2)i^(X + e/2)) . (73) 

^nuci fiu the full ( — h) self-energy for the nuclear pro- 
cesses; the current j stands here for charged or neutral 
nucleon currents defined in (j70p . Quantum states and 
operators are taken in the Heisenberg picture. The sum 
in Eq. (1721) runs over lepton spins. 

In the graphical form, the general expression for the 
probability of the lepton (electron, muon, neutrino) and 
anti-neutrino production is depicted as 



2 nucl,/j,z^ / J L J y 




(74) 



The hatched block has the meaning of the ( — h) self- 
energy of the virtual Z or W bosons which convert in a 
lepton / and an anti-neutrino v. The block determines 
the gain term in the generalized kinetic equation for the 
full (— -I-) Green's functions of the virtual boson. If one 
integrates over e//i (closes the e//i line in diagrams), one 
recovers the gain term for the i> in the charged current 
processes, see Refs. [47l. fl07L fl2^. This circumstance al- 
lows to use the given method in different non-equilibrium 
problems, like in description of the neutrino/ antineutrino 
transport. Note that the generalized kinetic equation for 
V is greatly simplified, if conditions for the qu asi-particle 
approximation are fulfilled, see [isl . ll07L Il25^. 

The integration over the lepton phase-space can be per- 
formed separately from the calculation of S,7u^i If we 
introduce the leptonic tensor 

TZM)- j mY.{^^n5^'Hq~qi~q.). (75) 



are to be constructed with the neutral and charged nu- 
cleon currents from Eq. ([7D)) . respectively. 

As we have discussed in the Introduction, the stan- 
dard calculation of the reaction rates is done with the 
help of summation of the squared matrix elements of re- 
actions, see [qO]. This is fully correct procedure, if one 
treats the processes perturbatively, i.e. provided there 
is a small expansion parameter. One nucleon processes 
are related to perturbative E ^ diagrams with only one 
the nucleon Gq^ Green's function in expansion of (|74l) . 
two-nucleon processes are related to the diagrams with 
two nucleon G^^ Green functions, etc. However, this 
procedure fails when applied to strongly interacting sys- 
tems. The description of even a one-nucleon process in- 
cludes infinite number of perturbative diagrams with the 
iViV-interactions, since the coupling is not small. Never- 
theless, one is able to separate processes using the quasi- 
particle approximation provided excitation energies are 
sufficiently low a; ^ ep (when the fermion width is 
small). Then, diagrams with one quasi-particle nucleon 
G ^ /F ^ Green's function describe the one-nucleon re- 
actions, with two G ^ /F ^ Green's functions describe 
two-nucleon reactions, etc. The calculations of the re- 
action rates based on application of the optic theorem 
formalism and calculations using the ordinary formal- 
ism of computing squared reaction-matrix elements yield 
the same results [ssl l46j . In general case, when parti- 
cle widths are not small, situation becomes much more 
involved. Then calculations using the squared matrix el- 
ements become invalid, and the only possibility to calcu- 
late the emissivity from the piece of matter is to use the 
closed-diagram technique. Below we formulate a general 
method and then discuss the quasi-particle approxima- 
tion. 



B. Diagrammatic decomposition in terms of full 
( — h) Green functions 



the reaction probability takes the form 



dW^.*_ 



.G2 



^J^--^^n.i>,..(^,<z)T-(.). (76) 

For evaluation of the emissivity in Eq. (157)) we also need 
the following expression 

jyytot (^2 

d^^d#^-^^'^^-^-^-v(^>^) 

xU;^(g), (77) 

where I — , fi" and 

= I u;, Y^i^^n s^'Hq -m- g.) . (78) 

spin 

The tensors T^^^ and Ty^^ are calculated in Appendix IB] 
Note that the self-energies, in Eqs. ^ and ((77)) 



1. Fermions with finite width 

The hatched block in Eq. (fM)) is the sum of all closed 
diagrams written in terms of full Green's functions. Ex- 
ternal ( — h) signs mean that each diagram in the se- 
ries contains at least one ( — h) nucleon Green's function 
(G and additionally F ^, for a system with pairing). 
The latter function is especially important. It obeys the 
Kadanoff-Baym kinetic equation. Various contributions 
from {X} can be classified according to the number N of 
exact G ^ /F ^ nucleon Green's functions (lines in the 
diagram) : 



I V 



(79) 
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The quasi-particle approximation for fermions can be uti- 
lized only if the energy of radiating quanta w is larger 
than the nucleon width (a; ^ e* ^ T ^ ^nX i^" 
equality w ^ T ^ ep rnust hold, see Refs. [sH . |47| . In 
this case the contributions of various processes encoded 
in a closed diagram can be made visible by cutting the 
diagrams through the (H — ) and ( — h) lines. The cut 
means taking off the energy integral provided the spec- 
tral functions of fermionic quasi-particles can be reduced 
to the (5-functions. This restricts the fermion energy to 
an in-medium mass shell. The N = 1 term describes the 
DU process, and iV = 2, the MMU and MNB processes. 

In general case, when the fermion width cannot be ne- 
glected, the cut through the ( — h), (H — ) lines has only a 
symbolic meaning. Nevertheless, the separation of the di- 
agrams according to the number of the full {—+) Green's 
functions proves to be helpful also in this case [47[ . Note 
that now each diagram in ([7^ represents a whole class 
of perturbative diagrams of any order in the interaction 
strength and in the number of loops. 

The full set of topologically distinct skeleton diagrams 
for S ^ written in terms of full ( — 1-) Green's functions 
can be explicitly presented as a series in TV [47|. For 
A'' = 1 and = 2 we have 



(80) 




For = 3 we have, for example, contributions of the 
type 



3)- 




H ^^^VWMi^^WW^f^^^ - 



(82) 



Note that for > 3 there appear multi-cut diagrams 
(see the last explicitly presented diagram in ([5^ ). The 
NN interaction block in Eqs. 



T8T|) and §^ is the full 
block containing the vertices of one particular sign, e.g.. 




(83) 



and the analogous equation for the (-I--I-) block. Since 
only the same-sign vertices are permitted in Eq. ()83|) . 
no (±,t) [i-e., (+,— ) and (—,+)] lines appear in these 
diagrams. The thick double-wavy lines stand here for 

the exact boson ( ) Green's functions or an iterated 

two-body potentials: 



(84) 



The full dot in the vertex is the renormalized in-medium 
vertex, which includes all diagrams with one sign, i.e. 
it is irreducible with respect to the full (±=f) nucleon 

Green's functions. This means that it contains only ( ) 

or (-t--l-) Green's functions. We denote such vertices as 




and 



e-g-, 




(85) 



Here we assume that the bare vertex is time-local, i.e. it 
carries on only one Keldysh index instead of three indices 
in general case. Note that full one-sign Green's functions 
entering Eqs. (|83l84l85p . would contain alternative sign 
diagrams, if they were expanded in perturbative series 
with respect to the bare Green's functions. To simplify 
discussion, we do not include in Eq. ([85)) the processes 
when the weak interaction (dashed line) is coupled di- 
rectly to the intermediate pion line due to the tttt — >■ Dl 
processes, see the first diagram in Eq. (|90p below. 

We did not show the direction of fermion lines in the 
diagrams, since it can be picked out at will in closed-loop 
diagrams. Once some direction is chosen, the arrows in 
the diagrams (|83l84l85p follow. 

For a theory of fermions interacting with bosons the 
contribution with the fewest number of external parti- 
cles is three (rather than four as for processes described 
with the Boltzmann kinetics). Indeed, the cut through 
the one-loop diagram in Eq. ()80|) shows that in dense 
matter an off-shell fermion can decay into a fermion plus 
a boson an vice versa. For these processes it is impor- 
tant that all particles have a finite width in dense matter. 
The formation and decay processes which are forbidden 
by the energy-momentum conservation in the free space, 
can occur in the dense matter without principal restric- 
tions. Therefore the most important term in the series 
([75]) is the first (one- loop) diagram (IM)) . which is posi- 
tively definite, and corresponds to the first term of the 
classical Langevin equation, for details see [47j . 

The classical Langevin process deals with the propaga- 
tion of a single charge (say a proton) in a neutral medium 
(e.g in the neutron matter). Therefore for this case only 
those diagrams occur, where both photon vertices attach 
to the same proton line. In the quasiclassical limit for 
fermions (with small fermion occupations nf <C 1; in case 
of equilibrium matter n/(po) = [exp((po ~ + 1]~^) 
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all the diagrams 



+ 




with arbitrary number of ( — h) lines are summed up lead- 
ing to the diffusion result, see Ref. [13] for detailed dis- 
cussion. Each of these diagrams with the n vertical in- 
sertions corresponds to the n-th term in the Langevin 
process, where hard scatterings occur at random with a 
constant mean collision rate F/ = — 251]^ . 

For particle propagation in an external field, e.g., for 
the scattering on infinitely heavy centers (proper Landau- 
Pomeranchuk-Migdal effect), only the one- loop diagram 
remains, where the fermion line is given by 



2. Quasiparticle approximation for fermions 

The one-loop diagram (IM)) calculated with the quasi- 
particle fermion propagators determines the one-nucleon 
reactions: the DU reactions and the PBF (in case of the 
superfluid matter) [H, [H, [gI] . The contribution to the 
DU process vanishes for n < . 

The two-nucleon processes are encoded in the N — 
2 term in Eq. (|79p and are given in the quasi-particle 
approximation by the diagrams 



+ 




since one deals then with a genuine one-body problem. 

In the general case of a non-equilibrium system all 
above equations, being presented in the Wigner repre- 
sentation, are very cumbersome. To simplify the prob- 
lem one may use the gradient expansion in the convolu- 
tions of two-point functions, see Eq. (|A11[) . In general 
case, the Wigner transformation will produce an infinite 
tower of nested gradient terms. Assuming that a piece of 
non-equilibrium matter under consideration evolves very 
slowly in time and space, one may keep only first-order 
gradient terms. First-order gradient terms in the expan- 
sion of the collision term Ccoi = S ^G^ — G ^ 
are attributed to the memory effects in the generalized 
kinetic equation |l07| . In the standard derivation of the 
kinetic Kadanoff-B aym equ ation one simplifying us ually 
drops these effects [13, Ill2j | . As has been shown in [l07j | 
the memory terms are of the same gradient order as 
other terms in the Kadanoff-Baym equation and should 
be kept, because of the local part of the collision term 
is also of the first gradient order (since Ccoi being zero 
in the thermal equilibrium state). However, in the given 
paper we are interested only in the calculation of the 
production rates in direct reactions from a piece of the 
non-equilibrium matter, which are fully determined by 
the quantity E ^. Since S ^ ^ in the thermal equi- 
librium, the gradient corrections to it are small and can 
be neglected provided the given non-equilibrium piece of 
matter evolves very slowly in time and space, that we 
further assume. Therefore, in further we will calculate 
only the local part of the E ^ term. 



Note that the first diagram in (|88)) is not allowed in terms 
of the full Green's functions with the width (compare 
with (|8T|) ) but it should be explicitly presented in the 
quasi-particle picture, see [35", After the cut over 

(— -|-), (-1 — ) lines the diagrams (|88| arc separated by two 
pieces and correspond to the processes 




(89) 



shown here for the paired potential interaction. 

In case when the TV TV interaction amplitude is mainly 
controlled by the soft pion exchange in the reaction chan- 
nel under consideration (for n > tiq), instead of one 
has 



p 




Ref. [Ill calculated the rate of MMU and MNB processes 
taking into account in-medium effects for the case of non- 
superfiuid nucleon matter. Evaluations [H, [H, have 
shown that the dominating contribution to MMU rate 
comes from the first two diagrams of the series ([50]) . 
whereas the third diagram, which naturally generalizes 
the corresponding MU (FOPE) contribution, gives only 
a small correction for n > ng. As is seen from comparison 
of Eqs. ((Ml) and ^U^, the first diagram dHHl) is absent. 
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if one approximates the nucleon-nucleon interaction by a 
two-body potential. 

The diagrams that can be cut into more than two 
pieces (e.g., see the last explicitly presented diagram 
in Eq. (1821) ) are proportional to powers of independent 
{L '")^+"(L^ )™ loops, m,n are positive integer num- 
bers, whereas the diagrams for a two-nucleon process 
have only two L ^ loops, and they decay after the cut 
into two parts. 

In Ref. [131 it was shown how for the processes with the 
radiation of soft quanta one can simply incorporate the 
effects of a finite fermion width into the results calculated 
within the quasi-particle approximation for fermions (i.e., 
the fermion width F/ = — 23Ey? is put to zero). For this 
purpose it is sufhcient to multiply the quasi-particle re- 
sult by a pre-factor. For example, comparing one-loop 
result at a non-zero value of the nucleon width F^ with 
the first non-zero diagram in the quasi-particle approxi- 
mation (when one puts F^ — >■ in the expressions for the 
Green's functions) one gets 



Coiio) 




(91) 



QP means here the quasi-particle result. At a small mo- 
mentum q of the radiated quantum the correction factor 
is equal to 



(92) 



which removes the singularity of the quasi-particle pro- 
duction rate for small uj. This factor complies with the 
replacement oj ^ uj + iTn in the retarded Green's func- 
tion. Correction factors for the higher order diagrams 
also can be derived. Here we quote corresponding results 
for the next lowest order diagrams: 





(93) 



(94) 



with Co{bj) from and 



^ N 



(Cj2 



(95) 



In case T <C ep for typical energy of the radiation lo ^ T 
one has w » Fjv, since Fjv cx T^/ep, and thus Cq — 
Ci ~ 1, see 47[. In general case the full radiation rate is 



obtained by summation of all diagrams in the series (1791) . 



C. Resummation of the two- fermion interaction 
out of equlibrium. Bosonization of the interaction 

In this section we work out resummation of the two- 
fermion interaction amplitude starting from a bare in- 
teraction, which is local in time but not necessarily lo- 
cal in space. This is the generalization of the procedure 
performed in Ref. [13] for the point-like interaction, lo- 
cal both in space and in time. We construct the com- 
pact expression for E ^ self-energy in equilibrium and 
non-equilibrium cases and illustrate how the decomposi- 
tion with respect to the number of (H — ) and ( — h) lines 
works. In order to simplify the consideration we first 
study the normal matter and then perform generaliza- 
tions to the superfluid matter. 

Gonsider the particle-hole channel with the full two- 
fermion interaction amplitude determined by 




(96) 

with some particle-hole irreducible bare time-local inter- 
action Qq. Without the first-order gradient terms in- 
cluded, the diagrams in Eq. (IM|) correspond to the fol- 
lowing expression in the Wigner representation 



dp 



fl4 



(27r)4i 



(97) 



X g^'{p',p"; q) G'\p" + |) G^'ip" - |) Q'^'{p",p-. q) 

where each element is to be understood as depending ad- 
ditionally on the Wigner's X variable. Since in the local 
approximation exploited here this variable will enter all 
expressions as a common parameter, we will not write it 
explicitly in the expressions below. In Eq. (j97p . p, p' and 
p" are the relative momenta in incoming, outgoing and 
intermediate channels, respectively; q is the exchanged 
momentum in the particle-hole channel. Note that in 
case under consideration the bare interaction is diagonal 
in the Schwinger-Keldysh space, i.e 



0, ^0 



^0. 



(98) 



We introduce the products of the Green's functions 

V\p- q) = -I &\p + ql2) G'^^p - q/2) , (99) 
which are related to the bare self-energies as 



— 



d/ 

W) 



Tl,{p-q)V^{p;q)Tl{p-q). (100) 



Since here the bare vertices are assumed to be time local, 
they carry only one Keldysh index, t]^ = r* 5ij5ik ■ For 
the vertex independent on the fermion momentum p the 
self-energy reads 



^0 



T^L^^q). V\q) 



dp4 



-V^ip;q), (101) 
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where we introduced the loop- functions L*-' . 
If we formally extend the products ([M)) as 

l^Hp; q) ^ L^'ip,?'; q) = l^'{p; q){27T)U^^\p - p') , 

Eq. (|96|) can be presented as 

4 



S-(,',P;,)=Sf(p',ri+i;/(g;j-/^^ 



X p";<))i"'(p",p"';?)S„"(p"',p;«)^ (102) 



The integral equation (I102p can be interpreted as a ma - 
trix equation in the discretized momentum space |l26l |. 
The integration turns into a summation over the grid 
points with the appropriate weights. Then in terms of 
matrices Q, Qq and L, Eq. ([M|) takes the form 



(103) 



where dot products emphasize the matrix multiplication. 
For example, the bare self-energy reads in these notations 
as Eq-* = Tg • L*^ • Tq. Working with matrices we can 
proceed with the solution of Eq. (|103p which we rewrite 
as 



g— = go- + g— ■ l— ■ g^- + g-+ ■ 1+- ■ g^-, 
g++ = g++ + g++ ■ l++ ■ g++ + g+- ■ l-+ ■ g++, 
g+- = g+- ■ L— ■ g^- + g++ ■ l+- ■ g^-, 
g-+ = g-+ ■ L++ ■ g++ + g— ■ 1-+ ■ g++. (io4) 

Introducing the quantity, called in [47[ the residual inter- 
action, 



^res 



^0^^ • [1 



= [1 ~ g^^ ■ L^^]-' ■ gt^ , (105) 

we rewrite the above set of equations as 

g = g^cs + 
g++ = g+t + g^ 

g^- = g^ 

g-+ = g- 

Substituting g^ , g ^ from the last two equations into 
the first two equations and using the notations 

I -1 



L+- 


■ Grcs 


L-+ 


"res 


L+- 




L-+ 


•^++ 

"res 



Z?^ = [1 - LT± • g;. 



±± 



L 



g^ 



zr = [i-g?,f-L^^-gtt-L^^r' 

we arrive at the formal solution of Eq. (|103[) : 



(106) 



g— 




■Z-- 




^rcs 7 


(107a) 


g++ 


= 


■ zr 


= zt 


+ .C++ 


(107b) 


g+- 


= ^++ 

^rcs 


■ zr 


■L+- 


• g — 






= g++ 


■L+- 


■zr 


^rcs ■> 


(107c) 


g-+ 


= g — 


■Z-- 


■L-+ 


■g++ 






= S^cs 


■L-+ 


■zr 




(107d) 



This solution for the bosonized interaction describes 
propagation of effective boson, such as phonon, plasmon 
etc., in non-equilibrium systems. These effective bosons 
can be treated on the same footing as all other effective 
quanta. The phase space distribution of such bosons is 
given by the g'^ and g ^ Wigner densities. 

Eq. ([5]) for the retarded Green's function decouples 
from the other Eqs. ([7]). Let us demonstrate that the 
same occurs for the resummed interaction (|107p . We de- 
fine the quantity 



-L- 



(108) 



which, being integrated over 4-momentum p gives the 
retarded loop-function = L — L as follows from 
^6p. Similarly we define 



(109) 



Using relations (jA12|) we are able to prove that these 

quantities are connected by the relation [L^] ^ = like 
the retarded and advanced Green's functions and self- 
energies. For the case of an energy-independent bare 
interaction go, Eqs. (jA13IA14IA15l) imply the useful re- 
lation 



go ■ {L— - L-+) -go^go- - ■ go . 



(110) 



With the help of this relation we obtain from Eqs. (|104p 
and Eq. 



g— + g-+ = ^0 + [g— + g-+) ■ l^" ■ go, 

-g+- - g++ = go- {g+- + g++) • • , 

where we used that {g -L^^ —g ^ •L++ + t/ ^ -L^ — 
g • L^^) ■ go (g + g-^) ■ -go and analogously 

{g++ ■ L++ + g+-^L-+ + g+- ■ L— + g++ ■ Z+-) • go = 

— (S^"*" + g~^ ) ■ ■ go ■ Thus, we can introduce the 
retarded interaction amplitude 

= g— + g-+ = -g+- -g++ (in) 
= ^0 • [1 - L" ■ go]-' = [1 - ^0 • L'"]-' ■ go , 

expressed only through the quantity , which convo- 
lution with go, Hke go ■ ■ go possesses the retarded 
properties. 

D. Physical meaning of multi-piece diagrams 

In general case the total radiation rate is obtained by 
summation of all diagrams in ([7^ . Some of the diagrams 
shown, e.g., in the second line in Eq. (|82]) give more than 
two pieces, if being cut. Therefore, they do not reduce 
to the Feynman amplitudes. The role of these diagrams 
will be illustrated in the given sub-section. 
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'^RPA 



1. Non-equilibrium systems 
Consider the RPA-like set of the self-energy diagrams 

(112) 

Note that this is only a RPA subset of all possible self- 
energy diagrams. In the Wigner representation with the 
omitted gradient terms Eq. (|112p reads 

S^PA = ^0 • ■r^+T.^o- L'" ■ Q^' ■ l'^ ■ rl (113) 

kl 

According to Eq. ()112p the quantity Sj^pA? which deter- 
mines the production probability, includes now the fol- 
lowing terms 



Taking into account Eq. (|llll) we obtain 



yR 



^RPA 



■"RPA 



^ '0 



(117) 



Equations (fTTT|) and (fTTT)) express the retarded interac- 
tion and the retarded self-energy through the bare inter- 
action Qo and the quantity determined by Eq. (jlOSp . 

It is possible to present expression for S^p^ in a more 
compact form. Using relations 



z±± 


^rcs 


= [1- 




5 


L-+ 


■ g++ 




L+- 








= L-+ 




-+ 


L-+ 




= z-- 




1 



- -+ 



(118a) 

(118b) 
(118c) 



(114) 



we rewrite Eq. (|115p as 



^RPA ~ '''res ■ 



■ L' 



(119) 



We see that if one singles out infinite tower of terms with 
only one ( — h) loop, their summation will lead to the 
renormalization of left and right vertices in the one-loop 
diagram, see the = 1 term in Eq. ((80)) . In addition 
Eq. (|113p contains terms with many repeated (H — ) and 
( — h) loops, i.e. the multi-piece diagrams. However, the 
RPA series does not include many other terms with N > 
2, as it is seen from comparison with Eq. (|82|) . 
From Eq. (|113p wc write now 

y-+ ~ ' 

^RPA 



where renormalized vertices 



= To 


■ L-+ 










+ ^(T 


■ L— 


■Q— 


■ L 


-+ 




+ ^(T 


■L-+ 




■ U 


f+ 




+ ^(T 


■L-+ 


■g+- 


■U 


-+ 




+ ^(T 


■ L— 


■g-+ 


■ u 


f+ 


-0+ 



Using another self-ener, 



■"RPA 



(115) 



= ^(T 


■ 








+ ^o" 


■ L— 


■g— 


L— 




+ V 


■L-+ 


■g++ 


L+- 




+ ^0 


■L-+ 


■g+- 






+ 


■ L— 


■g-+ 


L+- 





(116) 



we can determine the retarded combination S^p^ = 
^RPA + ^RPA- The direct calculation yields 



■"RPA 



^RPA ^ 'O ^ 



[g- 



To-L^'-ig- 



+ g-+)-L' 
-g++ + g+- 



g-^yi^-r. 



'res ^ 'O 

= [1- 



[l-L 



±± 



Go 



±± 



-±±1 



(120) 



are the solutions of Eq. (|85|) (and of similar equation 
for the (-I-) vertex) with the omitted second term on the 
r.h.s. (witin the RPA approximation). With the help 
of Eq. (|119p expressed in terms of L*^ one can calculate 
the reaction rates associated with the processes described 
by Eq. (|112p (in case of non-equilibrium slowly evolving 
systems with small spatial gradients). 



2. Equilibrium systems 

In equilibrium the expressions derived above can be 
simplified considerably and expressed through the real 
and imaginary parts of the function L^: 

go ■ L-+ ■ go = 2inbiuj)^go ■ • ^o) , (121a) 
go ■ ■ go = 2i[l + nb{uj)]<^{go ■ Z« • ^o) , (121b) 

go ■ L ■ go = go ■ ■ go 

+ 2mfc(Lj)5(go • ■ go) , (121c) 

go ■ L++ = -go ■L^'-go 

+ 2iK(cj) + l]3(go • • ^o) , (121d) 



where nf,(w) 



o/T 



1] ^ are the equilibrium boson 



occupations. These relations are derived in Appendix [Cl 
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We note that the self-energy (|119p can be written as Eq. (|117p we obtain 



^RPA 



'0 1^0 



■[l~g^+-L++]-'-r+, (122) 

where with the help of the equilibrium relations (|12ip we 
express 

2inb 



i + ?;(2n,, + i)y 



r++ T+- r + V r 

•^"=l-z(2n, + l)y 

[1 - g+t ■ Z++]-i . go =[1 - * (2nfc + 1) Y]-' 

■ [i-g^- -^L"]-' - go 

(123) 

through the common matrix 

Y = [i- g^- ■ 3fiZ^]-i • g^- ■ 3Z« . 

Since functions of the matrix Y commute with each other, 
we can simplify Eq. (I122p as 

Sr+a ^2in, tq- -g^^ • r • [1 + y • 
■ [i~go-m."r^ -4 

- 2 ^ nf, • g^^ ■ [ga ■ + (1 - • 5RZ^) 

• [^0 • 5Z«]-i • (1 - go • SRi'')] • T+ . (124) 

It is known that in the equilibrium the production rate 
can be also calculated with the help of the retarded self- 
enegy. Note that in this case the expression (I117P for 
^RPA '^^^ be also obtained from the direct summation of 
the series of diagrams 



(125) 



where all quantities are defined for T 7^ within the 
standard Matsubara technique (using discrete frequen- 
cies ioJn)- After the common replacement iujn oj + iO 
one obtains the analytical continuation to the retarded 
function. Applying standard rules to take the real and 
imaginary parts of the inverse of a complex matrix^ from 



^ If for a complex matrix M = 5RAf + iSAf, there exist an inverse 
matrix then the real and imaginary parts of the inverse 

matrix are equal to 

1 



sR(M-i) = 5RMS(Af"M . 



3A/ + 5RM^5RM 



(126) 
(127) 



5EgpA =To • ^o"^ • 3[1 - go ■ L"]-' ■ To 

=ro ■ go' ■ [go ■ ^L" + (1 - go5RZ«) 

• [eo3Z«]-i • (1 - ^oSRZ^)] • to . (128) 



Comparing Eqs. (I124p and (|128p we recover the standard 
equilibrium relation 



(129) 



see (jC4p . The sign minus compared to the corresponding 
relation for L, appears here since the self-energy includes 



vertices r. 







To ■ 



Thus, we have shown that in case of equilibrium sys- 
tems the reaction rates can be found either by using Mat- 
subara technique and then recovering 3E^p^ given by 
diagrams (I125p . or equivalently by summing up series of 
diagrams (|112l) within the Schwinger-Keldysh formalism. 
We emphasize that in the latter case not only the N = 1 
term but also corresponding multi-piece diagrams must be 
included. 

Let us now analyze the contribution of genuine one- 
nucleon processes separated according to the G*^ count- 
ing introduced in Sec. lIVBl They are represented by the 
only diagram ([50)1 



JV=1 



(130) 



It produces the expression 



Sjvtl = ^rcs + • T+s = To • [^0 ] ^ • g^c 



■L-+[i-g+t-L++] 

Making use of Eq. (|123p we obtain 



-1 ^+ 



(131) 



:]^ti =2in6To- • g^^ ■ [(2nfc + 1)2^0 • 

+ ii-go-UL'')-[go-^L^]-' 
■ ii-go-^L'')]-' ■T+ . 



(132) 



We see that in contrast to Eq. ()124p there is an extra 
factor (2 n;, -I- 1)^ in the denominator of Eq. (|132p . Thus 
the relation (|129p holds for S^llj^ only approximately, 
e.g. for n{, < 1. 

This example teaches us that multi-piece diagrams 
may yield a contribution to the total production rate, 
beyond that is given by the purely one-nucleon dia- 
gram (I8OI) . 

Moreover, as we will show below, only for the rates 
calculated with Eqs. (fm]) . ([TT7| (and Eqs. (fm]) . (fT28l) . 
respectively) the condition of the vector current conser- 
vation is exactly fulfilled. It would be fulfilled only ap- 
proximately (e.g., for Uf, ^ 1) or even violated in general 
case, if the rates were calculated according to Eq. (jl31l) . 
Bearing in mind that in many cases it is important to 
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keep conservation laws on exact level, provided it is pos- 
sible, we may re-interpret which diagrams correspond to 
the one-nucleon, two-nucleon, and other processes: We 
will ascribe a diagram to the one-nucleon process, if af- 
ter the cut through the (±=f) lines it decays into two 
pieces with two fermion legs, supplemented by the cor- 
responding multi-piece terms. Diagrams producing two 
pieces with four fermion legs each plus the corresponding 
multi-piece terms, describe two-nucleon processes, etc. 
We stress that only taking multi-piece diagrams into ac- 
count one recovers exact conservation laws (like the vec- 
tor current conservation) in sub-sets of diagrams respon- 
sible to one-nucleon, two-nucleon processes, etc. 



not be the same. The matrix of the L functions reads 



r 7^11 


Tl2 

^21 


^21 


L 


t21 
^12 


t22 
^11 


t21 
^11 


L 


^12 


Tl2 

^11 


^11 


L 


t21 
L ^22 


t22 
^21 


t21 
^21 


L 




(135) 



E. Extension to a superfluid system 

In the system with pairing we have to deal with 
the larger number of interaction amplitudes and loop- 
functions L. There are altogether 16 quantities corre- 
sponding to the possible direction of the arrows. The 
uniform description can be achieved in the so-called " ar- 
row space" introduced by A.J. Leggett in Ref. [16|, where 
each element is labeled according to the direction of the 
arrows attached to it. We use label 1 for the arrow to the 
left ("+" in the Leggett's notations) and label 2 for the 
arrow to the right ("— " in the Leggett's notations). For 
instance the particle-hole interaction Qq and Q entering 
Eq. (I103P will bring the indices (t/o)ii- The 16 elements 
of the interaction amplitudes and the loop can be now 
arranged in the 4x4 matrix according to the basis 



(133) 



The bare interaction matrix combines the interactions in 
particle-hole (5o)iii hole-particle {Go)22^ particle-particle 
(^^0)22: hole-hole (5o)ii channels arranged in the di- 
agonal matrix 



(.Go)f2 



(134) 



The indices and v run from 1 to 4. Note that interac- 
tion in particle-hole and particle-particle channels must 



In terms of the matrices Qq, Qq and L Eq. (|103p can be 
written as 

gj^ = g^ + EEg^-^^^-(go)'^ (136) 

jS kl 

in double-index notations, where the Latin indices run in 
the Schwinger-Keldysh space, j^i = +, — and the Greek 
indices run over the basis (jl33p in the arrow space, /z, z/ = 
1,2,3,4. Since all above derivations were performed for 
the matrix object, the generalization of the results (|107|) 
and ()111|) is now given in terms of the nested matrices 

and (unni). 

The equation for the self-energy (|113p takes the follow- 
ing form 



kl 



(137) 



where Tq is the array in the arrow space 
= (to'''^t^,0,0). The result (fTTOl) will be 
expressed now through the renormalized vertex 



^ros = (t^'-', T-', T^^)-', r^^^^) , obeying the equation 



l^.s-[l-Qr-L^T" = zt. (138) 
which has the diagramatic representation as in Eq. ()46p . 

F. Application to the point-like interactions 

To present the matrix relations derived above in a more 
transparent form, we consider a model case of a momen- 
tum independent bare interaction, Qq , and bare vertex, 
Tq . The matrix Eqs. (|107p turn into the algebraic ones 

with the replacement Go ■ L^-' ■ Go ^ Gq L^'' ■ The expres- 
sions derived in this section become very compact in this 
case. 
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1. Non-equilibrium systems 

The series of diagrams (|125p for the is easily 
summed up to the known result 



yR 
^0 



^RPA 



1 - eoS« 



RPA 



(139) 
(140) 



In order to simplify subsequent diagrammatic represen- 
tation we introduced the new object G = G/tq. 
On the other hand, from (I107P we find 



G-+ = 



[1 - Go^r][l + eoS++] + Gl^t'^o^ 



(141) 



It is convenient to re-group diagrams for Sj^p^ in Eq. 
([TT2|) as 



(142) 



ooo 



Thus, Ti + satisfies the equation 

y-+ _ _ 



[l-C?„-So-][l-^o++S++] 
Using (|107p we find 



^RPA 



[1 - goSo"-][l + ^oSri + ^oS^-Sq 



(143) 



(144) 



Similarly, for Sj^p^ diagrams in Eq. (jll2p can be re-grouped as 



Thus 



RPA 



^RPA 



y- + 



1 — C/q So 1 — ^/q ^0 

[1 - ^JoSq— ][1 + (?oS++] + ao'S+-E^- 



1 - C/Q 



(145) 



Obviously S^p^ is found from ()145p after the simulta- 
neous replacements ( ) ^ (++) and ( — h) o (H — ). 

Also Spp^ ()145p and Sj^p^ can be recovered from p43p 
and (|139p with the help of the relations (jASp . 



^. Equilibrium systems 

Let us now consider an equilibrium system. Using the 
equilibrium relations between TIq and the retarded loop 
self-energies: 



y- + 
^0 



2i3E^(l + nb). 



■i3I]^<(2nfc + l), 
+ i3E^(2nfc + l) 



(146) 



and comparing (|140p and (|144p we see that they satisfy 
exactly the standard equilibrium self-consistency relation 
(HHl). Adding and ([Ti5]) we recover p^ . 

Now consider = 1 term 



[1-goEo— ][l + goS++]' 



(147) 



where we took into account opposite signs of the bare in- 
teraction with the -|- and — signs. Using the equilibrium 
relations (|146p we obtain 



y- + 



^1 - go3fiS«) 



-2iQ?I], ^ 



(^0)2^0+^0+" 



^1 - ^o^fiE^)' + [(2nb + 1)^0^2^]' 



(148) 



Replacing (|140p and (|148p in ()129p we see that the lat- 
ter condition may hold for A = 1 term only approxi- 
mately, e.g. for rifc ^ 1. Also E^ll-^ ~ Sj^p^ provided 

|(2nfc -t- l)^o3E^| < |1 - ^oSRE^I that holds for T < ep 
in most interesting for us cases. However as we have 
stressed, it might be of principal importance to keep ex- 
act conservation laws, like the condition of the conserva- 
tion of the vector current. As we show below, the lat- 
ter condition is exactly satisfied provided one uses Erpa 
rather than Y,isi=i. 
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V. RENORMALIZATION OF INTERACTION 

The results of the previous section show that the in- 
teraction Qq enters the resuhing expression for the re- 
summed interaction [Eq. (|107|) ] and the self-energy 
11"+ [Eq. (|119p ] together with the L or Z++ functions 
in specific combinations given by the residual interaction 
Gtt [Eq- (UnSl)] and th e residua l vertex [Eq. P^ ]. 
The structure of Eqs. (|105|120p suggests that if one is 
able to split an a priori complicated function L" into 
some part, which we call "known", L^^^„^^, and a re- 
minder iJcm = ~^known' thcu the quantities GUs and 
T^^gg can be expressed through the "known" part iJtnown' 



±± 



±± 

cn 



±± 
cn 



±± 



r 



known J 

V 

known J 



^rcn 



and the renormalizcd interaction and vertex 



^rcn 



[1 



^0 



= [1-^0^ 



±± 



L 



rem J 

±±1-1 



±± 



(149) 



(150) 



this jump indicates the presence of a pole in the fermion 
Green's function. Thus, the full retarded Green's func- 
tion in the momentum representation is given by a sum 
of a pole term and a regular part Greg, 



e — e„ -I- i 7 



+ G«(e,p). (151) 



Here and below the energy e = po — ep is counted from the 

2 _ 2 

Fermi energy and the kinetic energy is ep = ^ 2 m*'" ■ 
residue of the pole term, a, quantifies the A.B. Migdal's 
jump in the momentum distribution 



Of. 



(152) 



where one can put T — which is correct up to higher 
order contribution 0{T^ /e^ 

The in-medium mass of the fermion is given by 



1 / 1 asf(e^ 



(153) 



By a cunning choice of the "known" part one can ac- 
count in ijjno^n for the most rapid variations with the 
energy and momentum in the interval of interest. Then 
the renormalized quantities (|150p will possess a weak 
energy-momentum dependence and can be cast in terms 
of phenomenological parameters adjusted to some em- 
pirical data. For nuclear physics such a renormalization 
program was conducted by A.B. Migdal in his seminal 
paper [5|. 



A. Fermi-liquid renormalization 

Simplifying consideration we focus below on the de- 
scription of equilibrium systems. Then we may deal with 
only one, e.g. retarded Green's function. The other 
Green's functions G"", G""*", G+~, and G++ are ex- 
pressed in equilibrium through the retarded Green's func- 
tion, see Appendix [Cl 

At low temperatures of our interest (T < A <C 
£F,nj£F,p) neutrons and protons are only slightly excited 
above their Fermi seas and all the processes occur in a 
narrow vicinity of eF,n and ep.p- In such a situation the 
Fermi-liquid approach seems to be the most efficient one. 
The basic assumption of the Fermi-liquid renormalization 
is that in a fermion system there is some mechanism of 
single-particle excitations. For normal systems at T = 
this manifests itself in a jump in the particle momentum 
distribution^. According to the A.B. Migdal's paper Q 



There is a special class of Fermi systems in which the jump is 
absent even in the normal st ate. They are called singular Fermi 
liquids or Non- Fermi liqiuds Il27ll . 



and the pole width is because of the coupling to the two- 
particle-hole mode. 



7 = -3^S,^(e,PF,Af)/e'' ^ const, 



(154) 



for e <C eF,N- At finite temperature we have to replace 
<^F,N ^ eF,JV + 0(T2/e|_jv) andpF,JV ^ V2 niN 

but these corrections are small for T < A and can be 
neglected. 

The typical averaged outgoing neutrino energy (w ^ 
T) is larger then the nucleon particle width Tis[{e ~ T) ^ 
7T^ ~ T'^/e-p.N- Therefore, one can neglect the width 
in the pole term in Eq. (jl5ip and work within the quasi- 
particle Green's function 



Gn,n.r>.{p) — G^.q.p. P ) — 



iO 



(155) 



Following 6|] only the pole part of G^ is relevant for 
descriptions of processes happening in a weakly excited 
Fermi system. The regular part can be absorbed by the 
renormalization of the particle-particle and particle-hole 
interactions on the Fermi surface. 

In the momentum representation the particle-particle 
interaction (j38p depends on spins and momenta of in- 
coming (a, q/2—p] b, q/2+p) and outgoing (c, q/2—p'; 
d, q/2 + p') particles [T^""(§ -I- p, f - p; f -t- p', f - 

P")\cd,ab ^ [^~~(P'^''9)]cd,ab' "^^^'"^ introduce the 
total momentum of two particles q and relative mo- 
menta in incoming and outgoing channels, p and p' . 
The particle-hole interaction (jSOp depends on spins and 
momenta of incoming (6, p + q/2] a, p — q/2) and out- 
going (d, p' + q/2-, c, p' — q/2) particles and holes, re- 
spectively, [U (p + f , p - § ; p' 



i dc.ab 
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[U — Graphical Eqs. (I35|), (|49]) can be 
written for causal functions: 

f dV 

Tpp {p,p'; q) = V ip,p'; q)+ J ^^-^ 

X V-ip,p"; q) G--ip'i) G--i-pl)f^--ip",p'; q), 

(156) 



^ph iP^P'yl) = U {p,P]q) + 



(27r)S 



X g) G--{p",) G--{p'L)f^-(p",p'- q) , 

(157) 

where p'^ — {e ± uj/2,p" ± g/2). The integration over 
p" involves energies far off the Fermi surface. One may 
renormalize interactions V in (jl56p and U in (jl57p such 
that integration will go over the region near the Fermi 
surface and one may use simple quasi-particle expressions 
for the Green's functions. Integral over internal momenta 
can be written now as 

where for T — Q the energy momentum integral is as 
follows 



d<i>o = P / ^ / de 



(158) 



and the brackets stand for averaging over the momentum 
direction n = p/\p\ 



47r 



(159) 



with p 



being the density of states at the Fermi 



surface. If the interaction V is not singular at g — 0, then 
for small q and \p \ ~ pf,n — \p'\ the amplitude Tpp is a 
function of cos^pp' = {nn'), n' = p' /\p'\, and Eq. (|156p 
can be written as 



+ ^pp(f«(n',n")fp-p-(n",n)),,, 



(160) 



The information about the bare interaction and the 
structure of Green's functions in regions far from the 
Fermi surface is encoded in the effective interactions 



T^{n',n)=V {p'f,Pf) + 



4^" 



d^P 



(2^)4 r^*' 

xV--{p'^,p'^)G--{p'')G--{~p'')rHn'',n). 

The quantity f (A ^ ^ ^ ep) does not appear in final 
expressions, except for the equation expressing the gap 
in terms of the phenomenological Landau-Migdal param- 
eter, see Eq. ()188p below. The ^pp function in the quasi- 
particle approximation for the Green's function (I15ip . 



e — ep + iO signe ' 



(161) 



renders 

Ap = J G- q.p. i+p) G- q.p. i~p) ei^~ ep) . 

The similar renormalization can be performed with 
Eq. (|157p . Here we have to note that even if U 
changes weakly for small q, the particle-hole scatter- 
ing amplitude is a sharp function of q because poles in 
G~~{p+) G~~{p^) approach each other producing (5(e): 

~<5(e) J de G,7_q.p.(p+) G,7_q.p.(p_) + 6(p, g), 

and B(j), q) is a smooth function of p and q. Then, 
Eq. (|157l) can be rewritten near the Fermi surface, \p \ ~ 
Pf.n ^ Ip'I, as 

fpj- 9) =r"(n',n) (162) 
+ {f-{n',n'')A-,-[{n''-q)f;,-{n\n-q))n^,. 

The renormalized interaction is determined by the equa- 
tion 



T^{n',n) = U—{p'^,PY) + 



4^" 



(27r)S 

X U—{p'^J'i)B{pv,q^G)f'^{n",n), 
where we assume \qp\/'m*^ « oj « ^ and 

B{pF,q 0) = lim lim [G'^"(pf+) G~"(pF-)] |, 
PF± = (±a;/2,pF n±q/2) . 
The particle-hole loop function in Eq. ()162p is defined as 

-^pir("; d$oG-q.p.(PF+) G,;-.p.(PF-)- (163) 



We emphasize that the amplitudes f ^ and f " are not 
calculable within a Fermi liquid approach and should be 
considered as phenomenological quantities. 

The renormalized amplitudes F^''^ can serve as a bare 
interaction 

in the analysis performed in the previous section. 

The renormalization procedure outlined in this section 
can be applied to a superfluid system without any modifi- 
cation, if the system is only slightly excited, T < A ^ ep. 
The regular quantities are subtracted at T — 0. The as- 
sumption of thermal equilibrium, which we used, allows 
to deal with the only one Green's function, here with 
G . If we deal with slightly non-equilibrium system 
for excitation energies e* ^ A <C ep, the same renormal- 
ization procedure should be performed for G ^. Actu- 
ally, it is more convenient to chose the spectral function 
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A = — as another independent quantity, rather than 
G '", and perform relormahzation of A (G = Af, 
where / is the distribution function satis fying the gen- 
erahzed kinetic equation, see Ref. (l07j |). After the 
renormahzation one may deal only with the quasiparticle 
Green's functions and with the renormalizcd interaction. 



B. Landau-Migdal parameters for the nuclear 
matter 

Because of the diagonal spin structure of the normal 
Green's functions, the spin structure of the amplitudes 
reflects the structure of bare interactions V and U : 

[^Hn',ri)]^^^^^^Ti{n',n){ta2)dc{t<y2)ab (164) 
[f ^)] dc.ab ^o{n\ n) {ao)dc {'Jo)ab (165) 

In the graphical form the quantity corresponds to the 
empty block in Eq. ([TU)) . In the absence of the spin- 
orbit coupling the interaction is invariant under indepen- 
dent rotation in spin and orbital spaces. In this case 
P?('^),a0 _ piti^) ja0_ Oppositely, in a isotropic system 
the interaction containing the spin-orbit coupling is in- 
variant under the combined rotations in both spaces and 
the spin structure of t he interaction is more involved. 
According to Ref. [l28j the general structure of the in- 
teraction in the system with one type of fermions can be 
written as 

r^{n',n)^r^<j'o<Jo+T^ (<?'<?) 

+ [3 {a'in-n')){a{n-n')) - {a'd)\ 
X (n~n'f 

-^V'^{a'{nY.n'\){a{ny.n'\) . (166) 

Each pre-factor here is a function of {n'n). The spin 
index assignment is the same as in Eq. (I165p . Matri- 
ces (Tj — {crj)ab with j — 0, . . . , 3, act on the incoming 
fermions, while matrices a'j = {(Jj)dc act on the outgoing 
fermions. The similar decomposition can be written also 
for r^(n', ft) with the only replacements aj ia2 Tj and 
armaria!,. 

The renormalizcd amplitudes, r^''» in Eq. (|166p . where 
a = 0,1,±,T, K , can be expanded in Legendre polyno- 
mials 



The structure with F!^ is the tensor interaction consid- 



some of the static properties of nucl ear and neutron mat- 
ter was found to be very small |129{ . However they could 
play an important role in the co nditi on for stability of the 
ground state of nuclear matter 'l30] . The spin-orbit and 
tensor terms in F^ are important for the description of 
3P2 pairing [79j . A non-local contribution of the pion ex- 
change to F^ and its effect on nucleon superfluidity were 
studied in Ref. 132]. For the description of the majority 
of nuclear phenomena including the 15*0 pairing it is suf- 
ficient to consider only the terms with Fq*-^-* and F^. In 
further we consider only these terms. 

Following Ref. [sj we introduce the dimensionless pa- 
rameters 



fi 



T 



a'^{no) p{no) ' 
F« 

a2(no) p(no) ' 



9i 



9i 



1-1 



a?{m)p{no) ' 
F« 



(168) 



where /'s and 5's are constants. We considered the sys- 
tem with the one type of particles, e.g. neutron matter. 
In general case of the nuclear matter of an arbitrary iso- 
topic composition, we need to know the nn, pp and np 
interaction amplitudes. These quantities can be param- 
eterized as 

+ 3iiJV2('^'*'^2)(*o-2i?), (169) 

where Ni,N2 = n,p. The calculation of the Landau- 
Migdal parameters is a formidable task and the results 
vary essentially, depending on a calculation scheme and 
a model for the bare nuclcon-nucleon interaction, see, 
e.g., Refs. [zl, ,133 ,136.] . Another possible path is to 
try to extract the Landau-Migdal parameters from the 
analysis of phenomena in atomic nuclei. Starting from 
Ref. Q one traditionally presents the nucleon- nucleon 
interaction amplitudes in the form, cf. [l37j |. 

C 

r = — — r [/ (To c^o + g (o- 0-) 



a2(no) 

+ f'^T'fa'„<To+g'" (r'f) (a '<?)], 



(170) 



(171) 

The constant C = l/p{no) is introduced as a dimen- 
sional normalization factor. One usually fixes its value 
as a2(no) C = 300 MeV-fm^, see Making use of the 
parameterizations (jl70|) and (jl71|) one implicitly assumes 
that the Fermi-liquid renormahzation preserves isospin 
symmetry of the strong interaction. Then, instead of six 
independent amplitudes f%^j^^ and g%^]^^ for nn, np and 
(167) pp channels one deals with four amplitudes f'^,g'^ and 



The others follow from the relations 



ered in Ref. [l29l - [T3ll| . The effect of the tensor terms on 



J nn J pp J 'J 1 J np J pn ~ J J i 



gnn = gpp = 5" + ff'", gnp = g%n = 5" - g'^- (i72) 
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This assumption can be justified only for the nucleon 
matter with a small isospin asymmetry. For strongly 
asymmetrical nuclear matter, like the neutron star mat- 
ter, the application of the relations (I172p is questionable 
and should not hold a priori. 

We note that for the particle-particle channel we use 
the spin parameterization (|169l) different from that in 
Eq. (|171|) . The two sets of parameters are related as 



/4 = (/^-/'^)-3(/-3'^), 
(/« + /'«) + (/+.9'«), 



-9^ 



(173) 



9'o- 
Jo 



The values of the zero-th and first Legendre harmonics 
of /j ffi 9' 9're extracted from analysis of many data 
on atomic nuclei. Unfortunately, there are essential un- 
certainties in numerical values of some of these param- 
eters. These uncertainties are, mainly, due to attempts 
to get the best fit to experimental data in each specific 
case, modifying parameterization of the residual part of 
the NN interaction. Numerical values of the parameters 
extracted in Ref. are ~ 0.25, ~ 1, ^ 0.5, 
1. Calculations in Ref. [138| give the values ~ 0, 
; 0.5-0.6, ~ 0.05 ± 0.1, 5^," ~ 1.1 ± 0.1. In 
Ref. 139] the value 9pp-o was fixed by the data on the 
two-neutrino double $ decays and the single (3 d ecays , 
as 9pp.Q — 1. This is in favor of the choice of Ref. |138| . 
First harmonics are related to the value of the 

effective nucleon mass. The values gpp.i = — 5p„;i — 
—0.11 are estimated from analysis of the decay energies 
and the Gamov- Teller strength distributions in neutron- 
rich short-lived nuclides (140.] . In Ref. [141| the values 
~ -(0.25-0.33) and g'^^ ~ -(0.3-0.5) are extracted, 
whereas Ref. [O^l calculated ~ -0.47 and g^ ~ +0.46 
using Cogny DSI force. 

Pairing gaps depend on the density and are very sen- 
sitive to values of parameters in particle-particle channel 
because of expone ntial dependence on the interaction am- 
plitude, see Refs. 142 1 and 0,185}. In application to the 
pairing in neutron stars it seems to be preferable to use 
the values of the Landau-Migdal parameters adjusted to 
reproduce the pairing gap s obtained in micorscopic cal- 
culations, like in Refs. (To-ISTI. lOTII . 



VI. EQUILIBRIUM SYSTEMS WITH PAIRING 

AT T 7^ 

In Section IIVI we have demonstrated that set of the 
diagrams for the one-nucleon process rate built up with 
the non-equlibrium Green's functions can be rewritten 
as the RPA series of the retarded self-energies. For the 
latter we may exploit the standard Fermi liquid approach. 
Thus further we follow the lines of 



A. Green functions and response of a system with 
pairing 

As argued in the previous section, for T ^ ep we 
can use the quasiparticle approximation for the normal 
Green's function, once studied process occurs near the 
Fermi surface. 

We assume that the renormalization procedure is prop- 
erly done. Thus, we may deal only with the pole parts of 
the Green's functions characterized by the effective mass 
TO* and the residue a. 

Then, neglecting (T/ep)^ corrections, for the retarded 
Green's function we write 

G^ip) = G^{p) ao , G^ip) = . (174) 

e — Ep + IV 

The Green's function of the hole is then given by 

G^^(p) = G^^(p) do , G'''«(p) = G^(-p). (175) 

Recall, here p = (e,p). We will use the approximation 
^p+q/2 ~ Cp + vq/2, where v is the nucleon velocity at 
the Fermi surface, v ^ vp n{l + 0(T^/e|)). Actually, the 
denominators of the Green's functions are \uj ± {£p+q/2 ~ 
ep-q/2)\ ^ ep- Moreover, the terms (x vq may vanish 
under the angular integrations. Taking this into account 
we estimate that the neglected terms are at most of the 
order of (A/ep)^ <C 1 compared to the remained terms. 
Such corrections are usually omitted in most calculations 
within the Fermi-liquid theory for superfluids. 

We need the gap function nearby the Fermi surface, 
hence, it is a function of the direction, n, of the relative 
momentum of pared fermions 



As we argued in Section IIIII the spin structure of the 
anomalous Green's functions repeats the spin structure 
of the gap functions 

F^'Hp) = {F^'\p)ao + Fi'\p)aya,, 

f'-'\p) = ia2{F^^\p)ao + Fi'\p)a). (176) 

Since, as we will see the gap is a sharp function of the 
temperature, we should retain this temperature depen- 
dence omitting it in other quantities. 

In momentum representation Gor'kov's equations (|31|) 
for the retarded quantities render 

G«(p) = G^(p) + G^(p)A(i)«(p,T)F(2)fl(p)^ 

^ G::«(p)A(2)«(p,r)G«(p). (177) 
Since G„ cx ctq one easily finds the solution 
-.R,^, _ ^o[Gr(p)]-^ 



G'\p) 



[G«(p)GJt«(p)]-i-HA2(p,r)/a2' 
F('^'^^ip) = A('^'^^ip,T)Gt^p)G^ip) 

_ A(i-2)«(p,r) 

- [G«(p)GJl«(p)]-i + A2(p,r)/a2- 
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We denote here 

A2(p,r) = -a2iTr{A(2)«(p,T)A(i)^(p,T)} 

= a'[Al{p,T) + Al{p,T)]. (178) 

Then using that 

[G„b) C^M]-' = [-(e + + 4]/a' , (179) 

we arrive at expHcit expressions for the quasi-particle re- 
tarded Green's functions in the presence of pairing 



(e + i0)2 - £2 



(e - -Bp + iO) {e + Ep + iO) ' 



F«(p) = 



(e + z0)2 - K2 



a UpVp 



{e-Ep + iO) {e + Ep + iO) 
where the Bogolyubov's factors are 



2 _ + 2 _ ^p ^p 

P 2E.„ ' P 2K ' 



(180) 



and the quasi-particle excitations possess the gapped 
spectrum 



El^el + A^p^T). 



(181) 



The energy-momentum integration at arbitrary tem- 
perature is defined as 

j &<^>Tf{e,ep) (184) 

/ + 00 
dep/(e,ep) for T = 
-OO 



-nr. 27ri 



n— — OO 



dep/(ie„,ep) for T 7^ 



The singlet-pairing gap is determined by the Fq term 
in the particle-particle interaction, and the gap equation 
reads 

A(n) = -(r«(n,n')Ao(A(rr'))A(n ')),!' , (185) 
Ao(A) = j A^rGaip) GHp) ~ Cp), (186) 

where Go{p) — l/(e„ — ep) is the Matsubara Green's func- 
tion for the Fermi system without pairing (A = 0). 

Note that the same value Aq can be introduced as An = 
Jd^T (G(p) G'^ip) + F{p) F{p)) - ep), cf. Rel 
For vanishing temperature the direct calculation gives 



ApjA) 
p 



+00 



de„ 



de 



27rz e-ep + zO e"^ - Ej + iQ 



de 



p 



= In 



In 



2i 
A 



2 _u ^2 
e> A. 



(187) 



After the retarded Green's functions are known other 
Green's functions can be expressed through them, see 
Appendix [C] In the rest of the paper we consider the 



singlet ISo pairing, hence A 



:i,2) 



and we will denote 



= aW ^ A (2) 

T = reads 



A = Aq ^ = Aq''''. So the causal Green's function at 



In case of the attractive interaction /q < the fermions 
are able to pair in ISq state and the relation between A, 
^ and /q is as follows. 



A = ^ exp 



p("o) 



(188) 



G-{p)^ 



a{e + ep) 



e^ — E^ +i Osgne ' 
-a A 
" e^^Ej + i Osgne 



(182) 



In further we will calculate the production rate using 
the equilibrium relation (|129p . Hence we have to calcu- 
late the retarded self-energy at finite temperature. We 
can use the Matsubara technique with 



ni \ a(ey^ + ep) 



-a A 



e2 - £;2 



g2 _ ^2 



, (183) 



In the renormalization procedure we used that ^ ^ ep. 
However, fixing the gap, one usually puts C = bearing 
in mind the weak logarithmic dependence. For T ^ 0, 



AoiA,T)^a^p 



dep(l — 2nf) 



(189) 



where Uf is the Fermi distribution. 

The diagrams for the self-energy are shown in Eq. (|57p . 
In momentum representation they produce the following 



where e„ = i 7r(2n -I- 1), n is the integer number running 
form —00 to 00. 



^ Definition of the value Aq is here the same as in Ref. [Tst and 
differs in sign from that used in Ref. |16| . llTll . 
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equation: 



tices on n, oj and q ), cf. Eq. ((55 



S — So + Si , 

So = ( y" d$T io l) (g+ G_ to(^, g) 
-F+ F_ g) + (G+ F_ - G_) foK <z) 



Si = ^ y d$T Q) [G+ G_ ii (n, q) 

+F+ F_ g) + (G+ F_ - G_) fi(r?, q) 



The left vertices here are the bare vertices (j45|) after the 
Fermi-hquid renormahzation, hke in Eq. ()150| . Their 
spin decomposition is f"(ri, g) = to{n,q)aQ + ti{n,q)a . 
Doing calculations in Matsubara technique, we use nota- 
tions: G+ = G{pq+uj,p + q/2) and G_ = G{po,p — q/2), 
and similarly for F± functions. For continues frequencies 
we use the symmetrical 4-vector notations, i.e. G± ~ 
G{p± q/2) and, analogously, F± — F{p± q/2). 

The full vertices, t(){n,q) and ti{n^q), are determined 
by the diagramatic Eqs. (|46)) . In Fig.|4]two of these equa- 
tions are written explcitly in the momentum representa- 
tion. These graphical equations were first introduced by 
A.I. Larkin and A.B. Migdal in Ref. [15i] . 



B. Larkin-Migdal equations 



The right vertices in Eq. (j57p are the full in-medium- 
dressed vertices, which are functions of the out-going fre- 
quency w and momentum q, and the nucleon velocity 
V c^i fi, fi = p/p . Their spin structure is 



f(n,q) = tQ{n,q)(TQ + &ti{n,q) , (190) 

T''{n,q) = to{~n,q) (To + a'^t{-n,q), 

T'-^\n,q) = {tlj^''{n,q)ao + at}^\n,q))ia2, 

T^'^\n,q) = ia2 {t''o\n,q)(To + dt^^\n,q)) . 



After opening the spin structure of the diagrams in Fig.|4] 
we arrive at the following set of equations for to, Iq, t'^^ 
and to (for brevity we omit the dependence of the ver- 



^0 — to 



thu 

^0 ''0 



d<S>Tr^ G+G-to-F+F^t^ 



G+ F^ t 



(1) 



F+ G_ t, 



(2) 



d$T n 



F+ G"_ t 



h ,(1) 







.(1) 



,(2) 



d$T 



G+ F^ to + F+ G'l t'^o 



(191a) 

G'l G'L t'^ - F+ F^ to 

(191b) 

(191c) 



G'i F_ tl?^ 



G+ Git 



h ,(1) 



d$T r« 



VI I VI if) 



F_ t 



(1) 



F+ G_ to + G\ F_ t^o 



(191d) 



The similar set of equations for 3-vector vertices ti, i/*. 



t'^'' and t\^' is written with the only differences that 
Fg'^ is replaced by F"'^ and in front of all terms with i'* 
the sign must be changed, cf. Eq. (|56l) . For the sake of 
convenience we introduce brief notations, e.g.. 



^(2) 



d^TGi F_ = G'] -F^. 



(192) 



The details of calculations of these products within the 
Matsubara technique are given in Ref. 95] . For instance 
the useful relations introduced in Ref. ^17'] 



F+-G'l_ = -G'l-F^, 



-F+-G- 



(193) 



are recovered for arbitrary temperature. From ()19ip we 
can immediately find relations between the vertices t^"^ 

and 4^^- Taking the sum of Eqs. (fT?Tc)l and (fT91d)) 
and making use of Eq. (jl93p we obtain the homogeneous 
equation for the sum t^^^ + 1''^\ 

4'^ +4"' = -( / d$T r« [g+ Gt - F+ F 



which implies 



4^^ + 4^^ = o 



(194) 



(195) 



for frequencies relevant for PBF processes, which we will 
consider below. The latter relation justifies the param- 
eterization of the full in-medium vertices used below in 
Eqs. (l202bl) - (f202dl) . The same relation is vahd for t'[^^ 

and t I vertices. 

In their original paper ^ A.I. Larkin and A.B. Migdal 
presented Eqs. (|19ip in somewhat different form. They 
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FIG. 4: Equations for the dressed vertices f and f'^' from Eq. (I46p in the momentum representation. 



noted that the vertices for the holes, t'', can be obtained 
from the particle vertices with the replacement n — >■ —n, 

toi"'^l) ^ to{~n,q) , t'l{n,q) ^ ti{-n,q) . (196) 

Therefore, one can introduce the operator tt, which per- 
forms this change of n in the vertex 

TTtoin,q) = t'${n,q) , T:ti{n,q)^t\{n,q). (197) 



Because of relation (|195p . Eqs. (|191c|191dp reduce to one 
equation for the vertex 

lo — -to - H ■ 

Analogously we introduce ti — —t^^ = i j^'* . Then four 
Eqs. (|19ip for scalar vertices '0' and four equations for 
3- vector vertices ti can be cast in terms of four equations 



io-to = (ro(iWio + Mto))^^,, 
h = - {tI {{N + Ao) + 0(7r) to 
r'^{L{-fr)ti+ Mti 



ti - t 



U = 



(198a) 
(198b) 
(198c) 
(198d) 



We shall call this set of equations the Larkin-Migdal equa- 
tions. In Ref. [l5l| these four equations are further re- 
duced to only two equations with the help of the operator 
P (defined by Eq. (31) in Ref. [IBl)) which includes addi- 
tionally the change of the sign between Eqs. (|198ap and 
([raSc)) and between Eqs. (ITMbl) and (ITgpt . Functions 
L, M , N, and O are defined as in Ref. Il5|: 

L(n,g;7T) = G+-G_ - tt , 

M{n,q) = G+-F^ - F+-G- , 
N{fi,q) = G+-Gt+F+-F^~Ao, 

0(n,g;7r) = -G+- F^ - F+^G^ti . (199) 



We emphasize that Eqs. (I198P are valid at arbitrary tem- 
perature. The temperature dependence is hidden in the 
convolutions of the Green's functions (I199P . In Ref. [lH| 
the latter ones were calculated explicitly only for T = 
using the method of Ref. psj . The extension to T ^ 
is done in Ref. |95|, see Appendix [P] 



VII. CURRENT-CURRENT CORRELATORS 
FOR THE PBF PROCESS n^n + vv 

Below we demonstrate how to apply the above results 
to calculate the neutrino emissivity in the PBF process 
n ^ n-\- vv from the superfluid neutron star interior. To 
perform this calculation 0, ^ we start with the dia- 
grammatic presentation of the current-current correlator 
([57]) for weak neutral currents, where enter the vertices 
shown in Fig. 21 calculated the previous section. In 
Eqs. (|571) . dSH]) we used the contour Green's functions, 
whereas now we will treat the same Eq. (j57P for the re- 
tarded self-energy. We use the Matsubara technique at 
T ^ and then perform appropriate analytic continua- 
tion to obtain the retarded self-energy. Thus we recover 
the RPA current-current correlator. 

The bare vertices generated by the weak currents (|70)) 
are equal to 



Tyin-, q) 

~Vfi 

T%in, q) 



ey ey ^ 

— , Tvi = — 

a ' a 

~gA{TXi<^lo-rlo<^l), (200b) 

-^uo (^A ^ 



Here ey and are effective charges of the vector and 
axial- vector currents. For the vector current ey = 1 and 
for the axial-vector current ~ 0.8 — 0.95, as it follows 
from studie s of the G amov- Teller transitions in nuclei, 
see Refs. giniliil and references therein. The corre- 
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spending vertices for holes are defined as 



where the contributions of vector and axial currents are 



Ty {n,q) = [Ty{-n,q)] 



(201a) 



= -9a{- T^i k - T^,o ^^l) ■ (201b) 
We used here exphcitly that T^fi{—p,q) = T^fi{p,q) and 

The structure of the full vertices reads 

T'v — gv {tv,oIo ~ T^v^l) o-Q , 

rv = gv{4olQ-T'^J)ao, (202a) 

= -5v (tv^o'o - , (202b) 
ta = -gA {ta,i^Io - ta,o^I ) , 
r'A = ~gA{r'X,,a'^lo-rl,a''l), (202c) 

^ gA{TAA^lo-TA.Q^l)icr2- (202d) 

We rewrite the retarded current-current correlator cor- 
responding to ([72]) as 



XV (q) = gliKok-r'^^J) (204a) 

X {loXv.o{n,q)-Xv.i{n,q)l)),-^^, 

XA{q) = 9A{{rlilo-rlol) (204b) 

X {loXA,i{n,q)-XA,o{n,q)l))~ 

with the scalar and vector response functions 

Xa,o(", q)^ J d$T(G+ G_ Ta,0 - F+ F_ T^Q 

+{G+F^-F+G-)Ta,o), (205a) 

XaA{n,q)^ J d$T(G+ G_fa.l+F+F-<1 

+(G+^_-F+G_)?a,o). (205b) 



cluJ'l" = Xiq) ^ Xv{q) + XA{q) , (203) 



Then we apply the Larkin-Migdal Eqs. (|198p for the 
case of the weak-current vertices (|201l202p . For the 
vector-current vertices we use Eqs. (|198all98bp and for 
the axial- vector-current vertices, Eqs. (I198cll98dl) . Then 
we separate the parts proportional to the scalar Iq and 
the vector I and obtain altogether eight equations for the 
vector and axial-vector current vertices. These sets of 
equations are cast in the following form [93 [. 



Tafi{ri,q)=T^^o{Ti,q) + {r^{n,n') [L{n\ q-,Va..o) TaA^^' ^ q) + M {n' , q)Ta.o{n' , q)]) , (206a) 
fa.o(n,g)- - {Ti{n,n')[(N{n\q) + Ao)7a,o{n\q) + 0{n\q;Va,Q)Ta,o{n',q)]),^,, (206b) 

Tas{n,q)=f^^i{n,q) + {r^{n,n')[L{n\q;Va,l)faAn\q)+M{n',q)TaAn',q)]).r,,, (206c) 
?aAn,q)^ - {ri{n,n')[{N{n\q)+Ao)?aA{n\q) + 0{n',q;raA)TaAn',q)])n'^ (206d) 

I 



where a = V for the vector and a — A for the axial- vector 
weak currents. The notation for the effective interaction 
is Ty^ = Fq'^ and F"^'^ = F^^'^. Operators 'Pa,i are de- 
fined as follows 

Va,^ = {-lTPa,^^T, Z = 0, 1 (207) 

with parameters 

Pv.o = 1 , Pv.i = -1 , Pa,o = -1 , Pa,i = 1, (208) 

which are the eigenvalues of the operators Va.i, acting on 
the bare vertices 

Pa.O<0 = ^a.O<o, ^a.l = ^",1 T^a • (209) 

To proceed, let us for simplicity assume that F^ and F| 
contain only zero-th Legendre harmonics. From (j206bp 



we find 

, (0(n,q:Pao))^ 

where 

, ^ ri{N{n,q)U .211) 

l + ri[Ao + {N{n,q))n]' 

For the channel a — V, for which the gap equation 
1 = —AoVy is valid, we obtain r/y — 1. For the other 
channel, J?^ 7^ 1. Substituting Eq. (piUl) in Eq. (I^D5a|) 
we obtain 

ra,o{q) = la{q; Pafi) Tafi > 

j-Hq;P)^l-K{^a{n,q;P))n, (212) 
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where we introduced the notation 

£a{n,q;P) = L{n,q;P) 

{ 0{n,q-P))n 
{N{n,q))n 



M{n,q). (213) 



Solving the second pair of the Larkin-Migdal equations 
(I206cl206dl) we first note that for the constants and 
the angular averages on the right-hand sides of equa- 
tions do not depend on n. Therefore, the component of 
the bare vertex proportional to v is not renormalized in 
medium. However, in view of the identity 



{f{n,q)n)n = {f{n,q) {nq))fiq/q" 



(214) 



valid for an arbitrary scalar function f of n and q, the 
full vertices gain a component proportional to q. Thus, 
we decompose 3-vectors Ta,i{n,q) and Ta,i{n,q) into the 
parts proportional to the n and fiq — vectors and 
introduce new scalar form-factors 



TaAn.q) = 'r^,iri + T^J'l{q)ng, 



(215) 



with the bare vertex r^i = {n-r'^ i) . Action of the oper- 
ator Va,i on the vertices (|215p is given by 

Va,lTa,l{n,q) = PaSTa.in+ {- Pa,l) Uq T^'^l{q) , 

'Pa,iTa,i{n,q) = {-Pas)nqTl^l{q) . 
Then, from Eq. ((206d)) we recover 



{O in,q; -Pa,i))« 
{N{n,q))n 



e {0{fi,q;Pa,i){nnq))fi 

Ta . 



(216) 



{N{n,q))n 

From Eq. (|206cp . substituting there Eq. (|216p . we find 

ri'l = -fa{q;~PaA)r': 

X {Cairi,q;Pa,i){n ■nq))fiT^,i. (217) 
Here we introduced the quantity 

Ca{n,q;P) = L{n,q;P) 

AM{ri,q))n ,^ 
~ "^^ /ATt- W 0{n,q;P) (218 
{N{n,q))fi 

and used the identity 

{£a{n,q;Pa.i)),i = {Ca{n,q;Pa.i))n 



(219) 



which allows to use for the vector vertices the same func- 
tion 7a as for the scalar vertex . 

In terms of the loop-functions (|199l) the response func- 
tions (|205al205bl) can be expressed as 

XaA^^ q) = -^(^1 <?; 'Pa,o) Taflin, q) + M {u, q) Ta.n{n, q), 
Xa,i{n,q) = L{n,q;Va,i)Ta,i{n,q) + M{n,q)Ta,iin,q). 



Using solutions (I210p and (|212l) . for the scalar vertices 
we find 



Xa.o(^, q) = 7a(9; Pa.o) <o ^H, q- Pa,o) ■ (220) 

With the help of Eq. (|215p we construct 

Xa,i{n,q) = L{n,q;Pa,i)nT^.^ 

+ L{n,q;-Pa^i)nqT^'^l{n,q) 

+ M{n,q)nqTi'^lin,q). (221) 



Using solutions (|216p and (|217l) for the three-vector ver- 
tices we obtain 



XaA^^q) = {L{n,q]Pa,i)n 



_ , , {0{n,q;Pa.,){nnq)U ^ 
M[n,q)fj, (^N{n,q))n 

+ Ca{n, q; -Pa^) la{q\ -Pa,l) 

X {Ca{n, q; Pa,i){nnq))n nq t^_i , 



and, then, rewrite it as follows 



1 'a,l 



(222) 



XaA^^q) = 'r^,ila{q; -Pa.i) Ca{n,q; Pas) 

+ SxaA^i,q), (223) 
SXa,i{n,q) 

{N{n,q))ri 

+ la{q; ~Pa,l) <1 (Cairi, q; -PaA 

X {Ca{n,q;Pas)n)n 

- Ca{n, q; Pa,l) n {Ca{n, q; ~Pa,l))n 



Relations 



(t^xv.o -qxv,i)n = 0, 

9((gu) {ujxvfi -qxv,i))n = 



(224) 



proved in Ref. 195| ensure the transversality of the po- 
larization tensor for the weak vector current. Thus, we 
have demonstrated that the retarded self-energy given by 
the RPA set of diagrams ((58|) with the vertices shown by 
Fig. 21 complies with the vector current conservation. 

Below we exploit expressions for the following averages 
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{Xa,o{n,q))n = la{q; Pa,o) T^fi 

X (£(n,Q;Pa,o))n, (225a) 
{Xa,o(n,q) {qv))fi = Ja{q; Pa.,o) 

X {C{7i,q;Pa^o)iqv))n, (225b) 

{qXaA^,q))n = la{q\~Pa,l) 

X {Ca{ri,q;PaA){Tl^q))n. (225c) 
{vXaAn,q))n = {L{n,q;Pa^i) {v f"^^i))n 

- {Min, q) {V n,)U {nW;^, 

+ 7a(g;--Pa,l)rQ {Ca{n,q;~Pa,l) {vnq))fi 

X (£a(n,g;Pa,i)(<i^9))fl- (225d) 

In the region of the ISq neutron pairing in neutron 
stars occurs at low densities (< ng) and the nucleons are 
non-relativistic, y-p <C 1. Expanding C and L at small Dp, 
see Refs. [o^, Ha], we find that the correlation functions 
7a differ from unity only in the second order in wp l^h 
i.e. 



7a(g;P)«i + 0(r:^ a' pq' v'^/uj') 



(226) 



We obtain that in the expression induced by the vector 
currents both scalar and vector components, (|225a|) and 
(|225dp . contribute at the order v^, 



S(xy,o(", q)) 
^{vxv,i{n, q)) 



4g fp 
" 45^4 



2q^V^ 



eyap%T(0, 1:^,0) , 
^ eyop55T(0, uj, 0) . 



Working in the leading order in vp we have to put ~^ 1 
in view of Eq. (j226p . The function 177^ is introduced in 
Appendix [Pl From (jD6|) it follows that 



%t(0,w,0) 



uy/oj'^ - 4A2 



VIII. NEUTRINO EMISSIVITY IN THE 
NEUTRON PBF PROCESS 



Once the averages (|225p are known, expressions 
and (|76p solve the problem of the neutrino emission from 
the superfluid neutron matter via the PBF reactions. 
The emissivity is given by 



^2 roo puj 



(227) 



48^4 



/ d|g I 



uj q 



Ka{q): 



where Ka is the correlator function summed over the lep- 
ton spins and integrated over the lepton phase space 

Ka{q) = 1^/ d<I>i35]xr('Z)- (228) 



The integration over the lepton phase space can be per- 
formed analytically with the help of Eqs. (|B2p and (I224p . 
and we find, see Ref. [9J for details, 



Kv{q) 



-{q 

a 

X ^{Xv.fi{n,q) ~ vxv.i{n,q))n , 



(229) 



2 

KAiq) = ^■^^•<s[q^{vxA,iin,q))n 
a 

+ (3w^ - 2(f^) {xAAn,q))n 

- uj{qXA,i{n, q))n - ^ ((gw) XA.o(^, <?))«] ■ 

(230) 

For the neutron PBF emissivity on the vector current 
we obtain [94,] (for one neutrino flavor) 



sPBF ^ (On) 2 2^4 22 



1. 



Here 



.(On) _ 



/(z) 



15 7r3 



(231) 



(232) 



is the result derived previously in Ref. [35|, |68|, |9l| . 

Now let us turn to neutrino emissivity induced by the 
axial- vector current. In the expansion tip ^ 1 the leading 
term contributing to the emissivity is of the order Vp. 
Keeping only the leading terms we cast Eq. (|230p as 



2 2 2 -'2 



Ka 

x[l + (l-|S)-|] %T(0,a.,0). 



(233) 



The correlation factors 7a contribute at the sub-leading 
order ^ Vp, therefore we neglected these terms in the 
approximate expression ()233p . As the result, the neutron 
PBF emissivity induced by the axial- vector current is (for 
one neutrino flavor) given by 



,nPBF 



1 



11 
21 



2 2 2 (On) 
9 A ^A^F,n '^vv 



(234) 



spin 



The resulting emissivity is the sum of contributions (|23ip 
and 1(234)) . We stress that Eqs. ((23T|) and ((234)) are ap- 
proximate expressions obtained in the leading order in 
fp. General result looks more cumbersome but it is eas- 
ily recovered with the help of Eqs. ()225p . The latter 
equations are derived in [95] and here at arbitrary tem- 
perature. 



IX. CONCLUSION 

A.I. Larkin and A.B. Migdal extended the Landau's 
Fermi-liquid theory onto superfluid systems. In this 
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paper we re-formulated their approach for systems out 
of equilibrium. For that we used Schwinger-Kadanoff- 
Baym-Keldysh formalism. Important improvements of 
the Larkin-Migdal approach compared to the Nambu- 
Gorkov one are that the former approach allows to deal 
with strong interactions different in the particle-hole and 
particle-particle channels. These achievements have been 
used by A.J. Leggett who generalized the Larkin-Migdal 
approach to describe strongly interacting fermion super- 
fluids at finite temperatures and applied it to description 
of superfluid '^He. He used Matsubara diagram tech- 
nique. The use of the Schwinger-Keldysh diagram tech- 
nique allows to consider variety of non-equilibrium prob- 
lems. In application to nucleon systems, in general, the 
considered in this paper formalism can be applied to the 
paring in the states with an arbitrary angular momen- 
tum; it operates with various forms of a nucleon-nucleon 
interaction: scalar, spin-spin, spin-orbit and tensor in- 
teractions. As argued by A.B. Migdal the tensor forces 
mediated by the pion exchange should enhance (pion soft- 
ening) with increase of the nucleon density. Inclusion of 
this effect might be important in the case of the P-pairing 
in neutron star interiors. 

We considered the neutrino radiation from a finite 
piece of the nuclear matter bearing in mind the prob- 
lem of the neutron-star cooling. We used optical theo- 
rem formalism formulated in terms of closed diagrams 
with the full fermion and boson Green's functions and 
the full nucleon-nucleon interaction. The series of the 
diagrams is constructed with respect to the number, 
of the full G ^ fermion Green's functions. For simpli- 
fication we considered a system which evolves slowly in 
time and has small spatial gradients. This allowed us to 
perform the gradient expansion after the Wigner trans- 
formation and keep only gradient-independent terms in 
calculations of reaction rates. We demonstrated that in 
order to exactly satisfy the vector current conservation 
in the nucleon pair breaking and formation processes it 
is not sufficient to include only one = 1 term of the 
series, rather one needs to re-sum the RPA series includ- 
ing multi-piece diagrams. (The multi-piece diagrams de- 
cay in more than two pieces, being cut through ( — h), 
(H — ) lines). This demonstration shows, how one should 
separate one-nucleon, two-nucleon, etc. processes, in ac- 
cordance with exact conservation law of the vector cur- 
rent. Comparison of the RPA S self-energy with the 
A'^ = 1 contribution shows the accuracy with which one 
may deal, using only one A^ = 1 diagram. 

Then we demonstrated how the developed formalism 
allows to calculate neutrino emissivity from the piece of a 
warm nucleon matter in presence of the nucleon pairing. 
As simplest example we calculated neutrino emissivity 
in the neutron pair breaking and formation processes. 
These processes are of one-nucleon origin. To simplify 
consideration, we focused on the case of the ordinary 
ISq pairing of neutrons. More difficult is to calculate 
the emissivity of the two-nucleon (A^ = 2) processes, 
and N > 5 processes in the presence of pairing. The 



existing nowadays results for the reaction rates in nu- 
cleon systems with pairing are based on the so-called R 
phase-space suppression factors used to reduce the pro- 
duction rates calculated without pairing. Such an ap- 
proach can be used only for rough estimations. The for- 
malism formulated in the present paper is fully suited 
to properly perform the calculations. It is also inter- 
esting to search for new processes which might be open 
in the non-equilibrium and equilibrium medium because 
of the interaction between different reaction channels. 
These questions require a separate consideration. In the 
present paper we focused on the neutrino radiation prob- 
lem. However, the white body radiation of other quanta 
can be considered in similar way. 

The calculated rate S ^ can be considered as the gain 
term in the generalized kinetic equation for the virtual 
W/Z boson or for the anti-neutrino. For consistency then 
one needs to include first-order gradient memory terms 
into the collision term. 

Another important question is how to go beyond the 
quasi-particle approximation for fermions in strongly in- 
teracting fermion systems with pairing. 

Within the quasi-particle approximation for fermions, 
the formalism based on the Fermi-liqiud renormalization 
is developed. To quantify the results it remains to know 
the Landau-Migdal parameters. For the problems under 
consideration one needs to know them as functions of the 
density, isospin composition, frequency and momentum. 
The information extracted from analysis of atomic nu- 
cleus experiments is definitely insufficient for these pur- 
poses. Existing calculations of the Landau-Migdal pa- 
rameters are still incomplete. We hope that the present 
study will motivate further attempts to extract these pa- 
rameters. In spite of all difficulties, it seems to be the 
cheapest way to achieve understanding of many new in- 
teresting phenomena occurring in the strongly interact- 
ing fermion systems in the presence of the pairing. 

The direction of the research was shown in the works 
done in 50th-70th years of the XXth century. The pio- 
neering contribution to the development of the methods 
of the quantum many-body theory including the problem 
of fermion pairing belongs to Arkady Migdal. 
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Appendix A: Matrix notation 

The Schwinger-Keldysh contour in Fig. [T] consists of 
two branches: time-ordered and anti-time ordered. For 
a given space-time coordinate x the contour coordinates 
take two values, a;"*" and x~ , depending on the branch 
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of the contour. The closed real-time contour integration 
can be written as 



dx'^ ■ ■ ■ — dx • • • — / dx^ 



(Al) 



to 



where only the time limits are explicitly given and the 
spatial integration d^a; is assumed. The folding of two 
two-point functions defined on the contour reads as 



H(x\/)= I dz^ax\z'')G{z^,y'). (A2) 



Any two-point function ( = (TcA{x)B{y)j being a 

function of two contour variables x and y can be viewed 
as a matrix, which elements are defined in dependence 
on their belonging to the branches of the contour 



C'{x,y)^C{x\y^), ^,Je {-,+}. 



(A3) 



For Green's functions this convention produces the fol- 
lowing matrices 

.^ikf \_ / iG {x,y) iG ^{x,y) 
iG [x,y) - iG++(x,2/) 

V($(x)$t(y)) (r$(x)$t(y)); ' ^ ) 

where T and T are the usual time and anti-time ordering 
operators and the upper sign is for fermions, the lower 
one is for bosons. Eq. (IA4[) implies the following rela- 
tions between non-equilibrium and usual retarded and 
advanced Green's functions 

G"ix,y)^G—{x,y)~G~+{x,y) 
= G+-{x,y)-G++ix,y) 
= e{xo-yo)[G+-{x,y)-G-+ix,y)], 

G^{x,y)^G—{x,y)-G+-{x,y) (A5) 
= G-+{x,y)~G++{x,y) 
^Q{y^~x^)[G-+{x,y)~G+-{x,y)\, 

where 0(a;o — ya) is the step function of the time differ- 
ence. The similar relations hold for the self-energies 

= -j:+-ix,y)-j:++ix,y) 

= 6(^0 -yo)[- S+- {x, y) + y)] , 

^^ix.,y)=^—{x,y) + i:+-ix,y) (A6) 
= -j:-+{x,y)-j:++{x,y) 
= 0(yo - XQ)[S+^{x,y) - J:~+{x,y)] . 



The difference in signs compared to (|A5|) is due to the 
fact that S includes vertices, and — and -I- vertices differ 
by the sign. 

In terms of matrices defined as in Eq. (|A3p the contour 
integration in the folding of two functions (jA2p turns into 



the usual space-time coordinate integration applied to 
the product of matrix-valued functions 



E 

!=+■- 



dzC\x,z)v'^G^'{z,y), (A7) 



where 77*^ is the diagonal matrix with the elements i] = 
—77++ = 1. It takes into account the extra minus sign of 
the anti-time-ordered branch of the contour. 

Applying hermitian conjugation operation to the 
Green's function definitions (jA4[) we obtain the follow- 
ing relations 



[G~~{x,y)]^ = -G++iy,x). 
[G'^ix^y)]^ ^G^{y,x), 



(A8) 



and the analogous ones for the self-energies. 

Instead of purely coordinate representation for two- 
point functions it can be more convenient to perform 
the Wigner transformation and operate in the mixed 
coordinate-momentum representation. For any two-point 
function y) one introduces the relative coordinate 
^ — x ~ y and the middle coordinate X = ^{x -\- y) and 
makes the Fourier transformation from four-space coor- 
dinate ^ to four-momentum p: 

C{X-p) = / d^e^P«C^(X + C/2,X-e/2) (A9) 



with i,j e {-I-,—}. The quantities G^ and G ^ 
are called the Wigner densities in the eight dimensional 
{X,p) phase-space. In the mixed Wigner representation 
Eqs. read 

[G±T(X;p)]t = -G±T(X;p), 
[G—{X;p)\' = -G++{X-p), 
[G^iX-.p]]^ ^G^{X-p). (AlO) 

In particular, these relations imply that the functions 
i G^^ are always real. 

The merit of the Wigner transformation is that it al- 
lows to set up an approximation scheme for treating the 
system being not to far out of equilibrium. For slightly 
inhomogeneous and slowly evolving systems, the degrees 
for freedom can be subdivided into rapid and slow ones. 
Then, the variable ^ relates to rapid and short-ranged 
microscopic processes and the variable X refers to slow 
and long-ranged collective motions. A gradient expan- 
sion with respect to slow degrees of freedom can be ap- 
plied since the Wigner transformation converts any con- 
volution of two-point functions into a product of the cor- 
responding Wigner functions plus first-order and higher- 
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order gradient terms 



d^e 



aX;p)G{X;p) + j 



dC do 



p'=p^X'=X 

dC do 



dpt" dXf, dXf' dp^ 



(All) 



Including local and only first-order gradient terms, 
one derives from the set of the Dyson equations the 
Kaddanoff-Baym generahzed kinetic equation describing 
slow evolution of slightly spatially inhomogeneous s ystem 
of particles having non-zero mass widths, see [l^ llOTl j. 
Although the gradient terms are assumed to be smaller 
than the local terms they should be kept not only in the 
Vlasov part of the generalized kinetic equation but also 
in th e collision term giving rise to the memory effects 
[lOTt . 

With the help of the relations (|A10|) we can prove the 
convenient relation between the functions i'-' introduced 
in Eq. (|99| : 



[L—ip;q)]^ = -L++{p;q). 



(A12) 



From Eq. (|A5P follows, inter alia, that among four quan- 
tities, e.g., and G**, not all are independent, be- 
cause of the relation* G++ + G—-G+--G-+ = 0. The 
similar holds for self-energies + S +E ^ = 
0, see. Eq. (jA6p . Let us check whether the similar com- 
pleteness relation holds also for the L functions. Apply- 
ing Eqs. (|A5P and (|A6|) recursively we find 



ones — and the integration contour can be close in the 
opposite semi-plane. Thus the relation 



dpo{L^ 



L— -L+- - L- 



= 



(A15) 



holds for L's but only after the po-hitegration. 

Appendix B: Lepton tensor 

Here we calculate the leptonic tensors entering the re- 
action probabilities (|76I77I) . We will assume that lepton 
1 is massive oji = \/iTif + |(fiP, whereas lepton 2 is mass- 
less UJ2 = I (72 I- The trace of lepton currents is given by 

^(Z^Zt)^2TrWi7^^^,.(l-75)} 



spin 



+ i eat,l3yqi q2, 

The integral over the lepton phase space 
d?qi A^q2 



2lui 2lu2 



< 92 '5^(91+ 92-9) 



can be calculated in the standard way. The result is as 
follows: 
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(l + 2^) {2q^q'^ + q'g^^)-3mlg^'^' 



1-^).H, 



(Bl) 



i{L++ + L— - L+- - L-+) = 

y^R \r^A I — ^ — h 



'p-q/2 



p-q/2\ 



+q/2 



'-'p+q/2 %+«/2j '^p-q/2 



'^p+q/2 



Gp+q/2Gp-q/2- (^13) 



Note that the last two terms vanish after the integration 
over po, 



The lepton tensor in Eq. (1751) is given by 
(1 + 2^) (,^,''-,^g-) 



47r 



^i-<)e{u)e{q^-ml). 



(B2) 



'^PoGp+q/2 Gp-q/2 ^ / '^PoGp+q/2 ^p-g/2 



0, 



(A14) 



since the poles of both integrated functions lie in one and 
the same semi-plane of the complex plane — 5po < for 
the retarded functions and 3po > for the advanced 



Actually in generalized kinetics only two real quantities are inde- 
pendent, e.g. iG ^ and A = — 2QG^. O ther quantities, , 
G , G^"*" are expressed through them [l07l |. In equilibrium 
only one quantity is independent. 



Setting mi = we obtain expression for the T-f^ptil) valid 
for neutral currents (for vD). 

For the process occurring on the charged current the 
reaction rates enters the integral: 



d^gi d^g2 
2uJi 2uj2 



9i 92 [1 - ■ w)]<5^'*H9i +92-9) 



where second term appeared due to the Pauli blocking 
for the charged leptons (electrons or muons) . Here u is a 
four-vector of a collective motion of the medium, in the 
co-moving frame u = (1,0,0,0), qi ■ u = 9fu^. Using 
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two four-vectors q and u we write the general structure 
of J^''' as 

+ D u"" + E g'"' , (B3) 
or alternatively as 

'-^ q ^ ^uu ^ ^uq Ju Jq 

+ C,^3[;fu+Cgg^^\ (B4) 
where we use tensors and four- vectors 

3^ ^q^- (u ■q)u^, ^ q' ~ {u ■ q) q^ , 

which satisfy relations 

f^^Uf.u^ = t'^.^qt^qi. ^q^ -{u- qf , 

^q,ti = 3 g , i^^^ = 3 , jqU^ = juQfi = , 

jqQf^ = j>t^ = q^-{u- qf , 

J^3u,^. = {u-q) [{u ■ qf - q^] . (B5) 
Coefficients in Eqs. (IB3p and (|B4I) are related as 

A = ^Cqq — {Cqu — Cuq){u ■ q) , 

B = Cquq^ + Cuq{u-qf, 
C = Cuqq^ + Cqu{u-qf, 
D = -Cuu- {Cqu+Cuq)q^ {u- q), 

E = Cg+Cuu+Cqqq^ . (B6) 

Using Ec^. (IB5P from Eq. (IB4I) we derive: 



4' =3 (g2 C,, + + 1] Q2 (^^^ ^ _^ 4 ^ 

Q.^Q.i' — Jqq — ^Cuu Q + Cg q , 

^{.iQv — Juq — Cuq Q H~ Cg ^ ^ 

(/^ Qf_iUi/ — t/gti — Cqu Q ~\~ Cg r2 , 
^l_L^i/ — Juu — Cqq Q Cg , 

Q^ = iu-qf^q\ n^iu-q). 



(B7) 



Explicit expressions for all Cy -coefficients look clumsy 
and we do not present them. However, they are easily 
written through the coefficient Cg, which equals to 

_ ^{Jqq + Q Juu) Q J^l ~ ^{Jqu + Juq) /-r!Ci\ 

2(2172 + Q2) ■ (B8) 

From the definition of the tensor J'^" we can easily find 
that the various convolutions with the vectors q and u 
and the trace can be expressed as 

J/^ = - m\) lo , Juu = nh-l2, 



Jqq = 5(9^ ~ "^1) ^0 , Juq = ^{q^ - ml) h 
Jqu = \ {q^ + mi)[^lIo - Ii] , 



where there appeared tree scalar integrals 
d^^i d^q2 , 



o o '-{u-qi)"n{u-q)S^^\qi+q2-q) 
T:6{q^ - ml)e{uj) 

X j d.ujin{uji) uj^yj uj^ — 

max{mi ,u;_ (g)} 



^±{q) = 



u} {q^ — m\) 



±k1 



Iq'^ + mj 
4g4 ' 



(BIO) 



n = 0, f,2. Finally, we can construct the lepton tensor 
for charged lepton currents 

TZt = 8{I^'' + r'^-g^''ip 

where 

8 ( J'^'^ + J"^^ - .g^'' JU) = 2 C,, + 2 C„„ tr 

+ [Cuq + Cqu) ijil fq + + 2 Cg g^'' 

- {Cqq 3 g2 + Cuu 3 + {Cuq + Cqu) QQ^ + Cg 4) g'^'' . 

For the tensor Tj^p^. we can use the same expression and 
the relations for the -coefficients with the only replace- 
ments — and Jab — > Jab, CL,b — q, u, where 

- - ml) ii , Juu = ni2-h, 

Jqq = 5(9** - nif) Ii , ,7„g = \{q^ - ml) I2 , 
Jqu^^{q^+ml)[nh~l2]. (Bll) 



Appendix C: Equilibrium Relations 

The equilibrium Kubo-Schwinger-Martin relations be- 
tween ( — h) and (-1 — ) fermion Green's functions and the 
boson self-energies are 



G-+{p) = TG+-{p)e-'/^, 
S-+(g)=±I]+-(g)e-^/^, 



(CI) 



see Ref. '149], where e = po — with the chemical po- 
tential related to the conserved charge. In the case con- 
sidered in the given paper /i ^ for fermions except 
for neutrinos/antineutrinos which freely escape from the 
piece of matter and /i = for bosons. All the Green's 
functions can be expressed through the retarded and ad- 
vanced Green's functions 



(B9) 



G—{p) = G''{p)±zn{s)A{p), 
G++{p) = -G^{p)±znf{e)A{p), 
C-+{p) = ±zn{e)A{p), 
G+-{p)^~z[lTn{e)]A{p), 



(C2a) 
(C2b) 
(C2c) 
(C2d) 
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where A = 2 Q^G"^ = — 2 is the spectral density and 

n(e) = l/(exp(e/r) ± 1) . (C3) 

Analogously for the self-energies we have 

E— (p) = S^'(p)±jn(e)r(p), (C4a) 

E++(rt = -E-^(p)±*n(e)r(p), (C4b) 

= T«n(e)r(p), (C4c) 

S+-(p)=*[lT"(e)]r(p), (C4d) 

where T = 2 -2 31]^ is the width. 

Now, using (IA12[) we derive the equilibrium relations 
among the products of the fermion Green's functions L*-' 
and the functions introduced as = L — L ^ and 
= L-+ - L++ = [L'^'W see Eqs (Il08ll09l) . Making 
use of Eq. (jCll) we immediately find 

iL-+lp] q) = ,„ G+" ,„ G+T ,„ G"+ ,„ e^T 

= iL+-{p-q)e-T . (C5) 



From Eq. piU)) follows 

= -2z3(go • L"^ ■ Qo) ■ (C6) 
Combining Eqs. (|C5|) and (jC6l) we obtain 

2^3(go•i'^•eo) 



^0 • L-+ ■ Go = 



(C7) 



and recover Eq. (jl21ap . Equation (jl21b[) follows then 
from Eqs. (|C7I) and (jC5|) . From the definition of the 
retarded function and Eq. (|110p we have 

Go ■ L ■ Gq = Go ■ ■ Go + Ga ■ L^^ ■ Go , 

Go ■ L++ ■Go = Qo- L+- ■Go-Go-L'' -Go, (C8) 



which together with Eqs. (IC7[) and (IC5[) translate into 
Eqs. ([mc)) and (fmd)) . 

For completeness, we list also the fermionic quasi- 
particle Green's functions for a system with paring in 
equilibrium 

G ^(e,p) — 2nianf{e) 

X [ul6{e - Ep) + vlS{e + Ep)]ao, 

_ /\(1,2)( \ 

X [,5(e + Sp)-<5(e-£;p)], 
G+-ie,p) = -27ria[l-?v(e)] 

X [u2^(e - Ep) + v^Sie + Ep)]ao, 

A(l,2)(p) 



F^^^^^+-{e,p) = 2TTia[l-nf{e)]- 



4 £'2 



X [J(e_Sp)-<5(e + £;p)]. 



Here the quasi-particle spectrum is given by E^ = + 
A^{p,T) with 

A'(p) = -a'^Tr{AW(p)A(2)(p)} 

and Cp — {p^ — pp)/(2mj^) . The Bogolyubov's factors 
are 



11 =: ^ — ?I — 

P 2Ep ' P 2Ep ■ 



(C9) 



For the Green's functions for holes 

G''±T(p)=a2G±T(-p)a2. 

Expressions for normal systems are obtained from here 
for Vp — )• 0. 

Appendix D: The loop functions 



At zero temperature the loop functions ()199p were cal- 
culated in Ref. 143] using the Feynman method for the 
integral of the Green's function products 



L(n,q-P) 
M{n, q) 
N{n, q) 

0(n,q-P) 



p 



ui — qv 2 



2 i^ + qv 

2 ^ ^2 - 



a p 



ci' p 



4A2 

ijj + qv uj — qv 



4A 



4A 



P 



(Dl) 



Here the variable P can be an operator as in Eq. (|199p 
or (±1) value as in Eq. (|210p . The universal function g 
is given by 



2 / d$, 



F+F^ 
a? p 



[e\-El] [el -Ely 



where /d$o is / d$T at T = 0, sec Eq. e± = 

e ± and E± = £'p±q/2 • Evaluation of the integral 
yields 



arcsinhV — 1 z tt ^(z^ — 1) 
z - 1 2z Vz^ - 1 

tj2 — ((fw)2 



(D2) 



4A2 



> 1 , U = Up n. 



In various limiting cases one obtains 
g{z) ~ 1 - zV2 for \z\ < 1, 



g{z) -i 



2\z\V-l-z^ 



ln(2z) 



for 



2|»»| 



for z^ ^ -1, 



for — z — > 



oo, 



(D3) 
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and for 
we have 



8A 



<\oj- 2A| < 2A 



O 



2|z|V^^-2A V 8A(u;-2A) 



A{uj - 2A) 



(D4) 



At finite temperatures the Feynman method does not 
work [l46| and the Matsubara technique can be used in- 
stead. In Ref. [9^ it was shown that the functions L, 
M, N, and O can be expressed through one universal 
temperature-dependent function 



gTin,iuJm,q) = 2 J d^x 

+?° +00 

^ J '^^p J2 



F+F^ 

2 

(2^ p 



= 2A^ 



(D5) 



T 



[(ze„ + lUJ^Y - El] [(*e„)2 - El 



Here e„ — (2n+ l)7rr and uj„i — 2rmrT. After the 
summation over n we obtain 



After the replacement ium W-i- = uj + iO we obtain the 
analytical continuation to the retarded gx function in the 
w-complex plain. The expressions for the loop functions 
read 



M(n,a;,q) = -a p ^—r — grin.uj.qj, 



2A 

,r/ ^ 9 — (qv)'^ 



0{n , uj,q; P) — ci' p 



4A2 

w + qv cj — qv 



L{n , uj,q\P) — a p 



4A 

X gT{n,uj,q) , 
V q 



9T{n,i^,q) , 
P 



UJ — V q 



4A 



(gTin, {vq),q) 



gT{n,uj,q)) - ^ grin ,i^,q] 



(D7) 



gT{n,iuJn,q) = A^ 



{E+-E^) {nf{E^)-nf{E+)) 



E+E- {iujm? - {E+ - E- f 

{E+ + E^) (1 - UFjE-) - nf{E+)) 
E+E- {lu^f ~ {E+ + E-f 



(D6) 



For T — Q, grin, {vq),q) = 1 and the old result (jDl[) is 
recovered. 
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